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CN . Abstract 



We introduce certain spherically symmetric singular Ricci solitons and study their 
stability under the Ricci flow from a dynamical PDE point of view. The solitons in 
question exist for all dimensions n + 1 > 3, and all have a point singularity where 
^^ ' the curvature blows up; their evolution under the Ricci flow is in sharp contrast to 

^N , the evolution of their smooth counterparts. In particular, the family of diffcomor- 

phisms associated with the Ricci flow "pushes away" from the singularity causing the 
P I , evolving soliton to open up immediately becoming an incomplete (but non-singular) 

■^^ • metric. In the second part of this paper we study the local-in time stability of this 

dynamical evolution, under spherically symmetric perturbations of the initial soliton 
metric. We prove a local well-posedness result for the Ricci flow near the singular 
2 . initial data, which in particular implies that the "opening up" of the singularity 

persists for the perturbations also. 
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5 The Linear step in the iteration: Proof of Proposition 14.11 

5.1 Proposition 14.11 follows from Theorem 15. II 

5.2 Plan of the proof of Theorem \5A[ A new iteration 

5.3 Apriori estimates for rj 

5.4 Energy estimates for ^: A Galerkin-type argument 

5.5 Improved regularity and energy estimates for ^ 

1 Introduction 

The question of defining solutions of geometric evolution equations with singular ini- 
tial data is an interesting challenge and has been studied in recent years for a variety 
of parabolic geometric PDE. For the Ricci flow, a number of solutions have been pro- 
posed in various settings. Simon [T8] obtained solutions for the Ricci flow for C" initial 
metrics that can be uniformly approximated by smooth metrics with bounded sectional 
curvature. Koch and Lamm [17j showed existence and uniqueness for the Ricci-DeTurck 
flow for initial data that are L°°-close to the Euclidean metric. Angenent, Caputo and 
Knopf [3] considered initial data of neck-pinch typeo They constructed a solution to the 
flow starting from this singular initial metric, for which the singularity is immediately 
smoothed out. This can be thought of as a (very weak) notion of surgery in that the 
method of proof relies on a gluing construction to show the existence of such a solution, 
but not uniqueness. Cabezas-Rivas and Wilking [5] have obtained solutions of the Ricci 
flow on open manifolds with nonnegative (and possibly unbounded) complex sectional 
curvature, using the Cheeger-Gromoll convex exhaustion of such manifolds. 

More results have been obtained in the Kahler case and in dimension 2, where the 
Ricci flow equation reduces to a scalar heat equation; we list a few examples: Chen, 
Tian and Zhang consider the Kahler-Ricci flow for initial data with C^'^ potentials 
and construct solutions to the Ricci flow which immediately smooth out. The argument 
is based on an approximation of the initial potential by smoother ones. Finally, more 
results have been obtained in dimension 2 (see [W] for a survey): Giesen and Topping |14] 
(building on earlier work by Topping [20]) have given a construction of Ricci flows on sur- 
faces starting from any (incomplete) initial metric whose curvature is unbounded; these 
solutions become instantaneously complete and are unique in the maximally stretched 
class that they introduce. More recently yet [15], they constructed examples of immortal 
solutions of the flow (on surfaces) which start out with a smooth initial metric, then the 
supremum of the Gauss curvature becomes inflnite for some finite amount of time before 
becoming finite again. 

This paper considers a special class of singular initial metrics and produces examples 
of Ricci fiow whose behavior is different from those listed above. Our initial metrics are 
close to certain singular gradient Ricci solitons that we introduce separately in the first 
part of this paper. The solitons exist in all dimensions n -|- 1 > 3. Our main result is 
that for small enough perturbations of the singular Ricci solitons, the Ricci flow admits 
a unique solution, up to some time T > 0, within a natural class of evolving metrics 
which stay close (as measured in a certain weighted Sobolev space) to the evolving Ricci 
solitons. In other words, we obtain a local well-posedness result for the Ricci flow for 
initial data with the same singularity profile as our Ricci solitons. 



^In particular these initial data can form in the evolution of a smooth spherically symmetric initial 
metric, as demonstrated in [T][5]. 



The solitons that we introduce (and in fact the perturbations that we consider) all 
have SO{n + l,IR)-symmetry. In particular, the soliton metric at the initial time t = 
can be written in the form: 

S-sol = dx^ +'^{xfgE.^, 

where x € (0, +00), for certain steady solitons (these we call half-complete, since they 

are complete towards x — )• +00) and x G (0, 5), 5 < +00 in the remaining casecl; here 

g^n denotes the canonical metric of the unit n-sphere. In all cases the function ^(x) 

is a positive smooth function and moreover ip{x) — )• as x — )• O"'', with leading order 

1 
behaviour ^ ~ x^. In particular, the (incomplete) metric above can be extended to 

a complete C^ (in fact C^) metric at x = 0, but the extended metric will not be of 
class C^ . We remark that (in the half-complete case) our (singular) solitons are complete 
Riemannian manifolds towards -|-oo, with an asymptotic profile there that matches the 
Bryant soliton. For the rest of this introduction we discuss only the half-complete case. 

Our first observation is that the evolution of the singular solitons themselves under 
the Ricci flow is in sharp contrast with the behavior of their smooth counterparts. As we 
note in Section [2.31 below, the evolution of the Ricci soliton g^ox is given by a 1-parameter 
family of radiate diffeomorphisms pt : (0, +00) x S" — )• (0, +00) x S", t > 0, where pQ = Id 
and such that the pullback g{t) = Pt {qsoi) solves the Ricci flow 

dtg = -2Ric{g), g{0) := gsoi- 

However, the map pt is not surjective in this case. In fact, for each t > 0, pi(0, 00) = 
{m{t), +00) where m{t) > is non-decreasing in t. In other words the flow pt pushes 
away from the singular point x = 0. Thus, for each t > {M,g{t)) can be extended to 
a smooth manifold with boundary, where the induced metric on the boundary is that 
of a round sphere of radius lim^_^o+ i^iPtix)) > 0. One can then visualize the evolving 
soliton metric g{t) backwards in time: Starting at time t = 1 it contains the portion of 
the original soliton corresponding to x > m{t), and its boundary at x = m{t) shrinks 
down, as i — >• 0+, to the singular metric ^soi- 

The perturbation problem that we consider is still within the spherically symmetric 
category. In particular, the initial metrics we consider are in the form 



which is such that if we let 



g = dx^ + w{x)gs^ 






then ^ is in a weighted ff ^-Sobolev space 



^eH\0,+oc) I 4r + 4^d^<oo 



X'" x^ 



for a large enough constant a which we determine below. In particular, this assumption 
forces the perturbation g to agree, at t = 0, with g^oi (asymptotically as x — )• 0^) to a very 



^So in the general case the sohtons can be extended to a smooth metric with boundary at x = 5 
^"Radial" here and furtherdown means that the diffeomorhpism, for each t>0, depends only on the 
parameter x £ (0, 00). 



high order, thus controlhng the asymptotics of the curvature tensor near the singularity. 
We then show that (subject to suitable Dirichlet boundary conditions at x = 0, t > 0) 
there exists a unique evolving metric g{t), t G [0,T], solving the Ricci flow equation, and 
which stays "close" (measured in a suitable weighted H^-space, the weight depending on 
both the spatial and the time variable) to the evolving soliton metric; in particular it 
exhibits the same "opening up" behavior of the initial singularity. The precise statement 
can be found in Theorem 13. 1[ 

It should be stressed at this point that our work here does not have direct bearing on 
the issue of "flowing through singularities" that form in finite time under the Ricci flow, 
(as studied, for example, in ^), at least for closed manifolds. Indeed, it is well known 
that for such manifolds the minimum of the scalar curvature is a non-decreasing function 
under the Ricci flow; however the scalar curvature of the solitons we consider (and of their 
perturbations) converges to — oo at the singular point (x = 0). Nonetheless, the issue of 
proving well-posedness for geometric evolution equations with singular initial data is an 
important challenge in many parts of this subject. In fact, from a purely PDE point of 
view it seems natural to first attempt this well-posedness in the neighborhood of solitons 
(or more generally stationary solutions of the flow, the stationarity being understood in a 
suitable sense depending on the context). It is hoped that the methods developed herein 
will prove useful in such pursuits in the future. 

1.1 Outline of the ideas 

Now, we briefly outline the sections of the paper and the challenges that each ad- 
dresses. In Section [2] we introduce the (singular) spherically symmetric Ricci solitons 
that we consider. The study of these solitons (all of which are gradient solitons, i.e., the 
one parameter family of diffeomorphisms associated to them is generated by the gradient 
of a function) follows the method presented in [Gj Chapter 1], originally developed by R. 
Bryant. In the class of spherically symmetric metrics, the gradient Ricci soliton equation 
reduces to a second order ODE system, which can be transformed into a more tractable 
first order system in parameters (W, X, Y) via a transformation that we review in ()2.6p . 
Knowledge of the variables W, X, Y in the parameter y allows us to recover the metric 
component ip and the gradient u of the potential function cj) of (12. Sp in the parameter 
X. In the case of steady solitons, the system ()2.8p in fact reduces to a 2 x 2 system; see 
§2.21 We provide a description of the trajectories in the X, F-plane that correspond to 
our singular solitons and compare them to the Bryant soliton. In particular, we show 
there exists a 1-parameter family of singular gradient steady Ricci solitons; they are all 
singular at x = with the leading order asymptotics 

/ \ ^- V^ ~ 1 

V'(x)~x^ a;(x) ~ , n>l 

X 

and they are complete towards x = -|-oo, with the same asymptotic profile as the Bryant 
soliton. 

For completeness, we also consider the case of spherically symmetric shrinking and 
expanding solitons; these correspond to trajectories of the 3x3 system (j2.8p . In this case 
we can again derive the existence of solutions of the above type, but only for x G (0, 6), for 
some small 6 > 0; the reason is that the behavior of the trajectories for large parameter 
y is not well understood. While the above solitons have been constructed over the 
manifolds M x S", it would perhaps be natural to seek similar examples in the more 
general cohomogeneity-1 category, studied by Dancer and Wang, [TOl [TT| [T2]. 



In Section [3] we introduce the perturbation problem we will be studying in the rest 
of the paper. We consider spherically symmetric initial metrics of the form 

For such initial data, the Ricci flow equation can be written (after a change of variables) in 
the equivalent form (I3.5p of a PDE coupled to an ODE. The evolving Ricci soliton metric 
(defined through a radial pullback) remains spherically symmetric and is represented by 
coordinate components x{xi't),ip{x,t), while the stipulated Ricci flow that we wish to 
solve for with the above type of (perturbed) initial data corresponds to two functions 
x{x-, t),ip{x, t). Since the singular nature of the initial data do not allow the system (13. 5p 
to be attacked directly, we introduce new variables which measure the closeness of x, V' 

to x,ip- 

More precisely, we define 

X ip 

Then the system reduces to (IS.lOj) . for which the Ricci soliton corresponds to the solution 
^ = 0,^ = 0. The coefficients of this system refer to the variable ip of the background 
evolving soliton, expressed with respect to its arc-length parameter s. What is critical 
here is that the coefficients are singular at {x,t) = (0,0); the precise nature of this 
singularity is essential in our further analysis. 

A first challenge appears at this point, which in fact is independent of the singularities 
of the coefficients. Indeed, it is related to the presence of the second order term ^ss on 
the RHS of the first equation in (|3.10p . Since the first equation is only of first order in 
C, this term would not make it possible to close the energy estimates for our system. We 
therefore introduce a new variable defined by 

_ (C + 1)^ . 

^ (^ + 1)2" 

The new system (j3.12p for r] and ^ involves only first derivatives of ^ in the evolution 
equation of rj and therefore can (in principle) be approached via energy estimates. It 
is unclear at present whether there is any geometric significance underlying this change 
of variables; it is in fact not a priori obvious that such a simplification of the system 
should have been possible via a change of variables. It is at this point that the spherical 
symmetry of both the background soliton and of the perturbations that we study is used 
in an essential way. 

Thus, matters are reduced to proving well-posedness of ()3.12p . in the appropriate 
spaces. We follow the usual approach of performing an iteratioiQ, by solving a sequence 
of linear equations for the unknows [r]'^'^^ , S,™"^^) in terms of the known functions (r/™", ^™') 
solved for in the previous step, and proving that the sequence {rj"^, ^™'), m G N converges 
to a solution (rj, ^) of our original system. 

We note that the usual approach would be to replace only the highest order terms in 
the RHSs of ()3.12p by the unknown function ^"^+^ and replace all the lower-order ones 
by the previously-solved- for rj"^ , ^"^ . However in the case at hand this approach would 
fail for any function space, due to the nature of the singularities in the coefficients. For 



^ In reality a contraction mapping argument, although we find it more convenient to phrase our proof 
in terms of the standard Picard iteration. 



example, as we will see the coefficient % in the potential terms contains a factor of -^, 
where s{x,t) is the arc-length parameter of the background evolving soliton. It turns out 
that the leading order in the asymptotic expansion of s^ near a; = 0, i = is of the form 

s^ ~x2 + {y/n- l)t. 

Consequently, the best L^ bound for -^ would be ^ < ^j which evidently fails to be 
integrable in time; this would result in an energy estimate of the form dt£ < j which 
cannot close. The remedy for this problem is to modify the iteration procedure according 
to (|4.ip . In this linear iteration the unknown functions ^"^+1^7^™-+! at the (m + l)-step 
also appear in certain lower-order terms associated to the most singular coefficients. 

Finally, we solve the system (j4.ip and prove that it defines a contraction mapping in 
certain (time-dependent) weighted Sobolev spaces H^{s) containing all functions 

T^ ^ + T^ ^-^ — rds < -hoo, 

(s2 + at)'' (s2 + at)°-i 

where we note that the weights depend on both the spatial and time variables x, t. We 
note here that the integration is over the entire half- line M+Ij but we use the length 
element ds which corresponds to the arc-length parameter of the background evolving 
Ricci soliton. In particular s{x,t) := pt{x); thus for all t > s{x,t) > 0,Vx > 0. 

The rather involved estimates in Section U] aim precisely to show that the parameters 
a and o" > can be chosen in a way to make the estimates close. The Picard iteration we 
follow is presented in ^4.H where we also assert that these linear systems are well-posed 
in the appropriate energy spaces, and can be solved up to a uniform time T > 0, with a 
uniform bound on the energies. Assuming that the well-posedness of the linear systems 
holds, the main claim is that these iterates form a Cauchy sequence. This is done in the 
subsequent subsections in 21 by first deriving bounds for the iterates in the weighted L^ 
spaces, and then in the appropriately weighted H^ spaces. 

In Section [5] we then prove the well-posedness of the linear systems ()4.ip : in fact we 
prove this for a more general type of linear equations. The approach here is to iteratively 
solve for a sequence of functions which solve a coupled PDE-ODE system (a heat equation 
in ^ and an ODE in r]). The ODE can be solved separately and the energy estimates 
follow readily. The heat equation in ^ is then solved by a Galerkin-type argument, proving 
that it admits an a priori weak solution (which a posteriori is a strong one) to which 
the energy estimates can be applied. This part is included for the sake of completeness, 
since coupled systems of this singular nature do not appear to have been treated in the 
literature. 
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^Recall that for the purposes of the introduction we are only considering the stabihty of the half- 
complete (steady) Ricci solitons. 



2 Singular spherically symmetric Ricci Solitons 

We will be considering metrics over M"+^ = (0, B) x S" (where B £ M-|_ oi B = +00), 
in the form 

g = dx^ +il)^{x)g^n, (2.1) 

where ip \s a, positive smooth function and g^n denotes the canonical metric on the unit 
sphere. 

Our first aim for this section is to obtain such metrics which satisfy the Ricci soliton 
equation and which are singular as a; — t- O"*". In particular we wish to construct a soliton 
metric which will extend continuously to x = with 'il:{x) — t- 0, as x — t- 0^, but will not 
close smoothly there (see ^ Lemma A. 2]). (In the case B < +00, we will consider ^(x) 
having a smooth limit at x = B). 

Recall the gradient Ricci soliton equation 

Ric{g)+W(p + Xg = A G M, (2.2) 

for a smooth potential function </> : M — )• M, depending only on x, (l){x,p) = (p{x). 
Generally, for metrics of the form ()2.ip , we can easily check that the Ricci tensor is given 
by 

Ric(g) = -n--dx^ + in - I - ipi) - in - l)ii>f)g^-a (2.3) 

and the Hessian of a radial function (j) by 

VV(/> = ipdx'^ + 'tpip<t>gnr^ , (2.4) 

where = ^- Therefore, setting oj = cj), equation (|2.2p reduces to a coupled ODE system 
of the form 



(2.5) 



^-0 + (n — l)'ip'^ — (n — 1) — iIjt^uj = Xip'^. 
Following [6l Chapter 1, §5.2], we introduce the transformation 

1 ^ ^ Jnin-l)W , , 

-uj + n^ 1" '/' 

along with a new independent variable y defined by 

dx , „, 

dy = -. (2.7) 



For the above set of variables, our ODE system (12.51) becomes 



d 
dy 



' W = W{X'^ -\W^) 
X' = X^-X + ^ + \{^-X)W^ (2.8) 

Y' = YiX^ - ^ - XW^) 

We readily check that for the steady (A = 0) and the shrinking (A < 0) cases, the 
equilibrium points of the above system are 



(0,0,0) (0,±1,0) (o,^,±Jl-i) 

\/n V n 



In the expanding (A > 0) case, we have two additional critical points (ib-i=, ^,0). 
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Remark 2.1. In the steady case, A = 0, the system Ii2. 8\) is in fact represented by the last 
two equations (W is a function ofX,Y). The well-known Bryant soliton corresponds to 

the unique solution of the system that emanates from the equilibrium point, {—t=, y 1 ~ :^) 
and converges to the equilibrium {X = 0,Y = 0) J^ Chapter 1, ^4]- 

In the present article we are concerned with the trajectories emanating from the 
equihbrium point (0, 1,0), for ah A G M (in our primary analysis). 
The linearization of (j2.8p at (0, 1,0) takes the diagonal form 

I I X-1 I (2.9) 

Note that for n > 1, all eigenvalues (diagonal entries) are positive, which implies that 
(0, 1, 0) is a source of the system. Whence, if a trajectory of (|2.8p is initially (y = 0) close 
to (0,1,0), i.e., 

\{W{0),X{0)-l,Y{0))\<e, 

for e > sufficiently small (indicated by the RHS of (j2.8p ). then standard ODE theory 
(e.g., see ^) yields the estimate 

\{W{y),X{y) - 1, Y{y))\ < ^/3ee'^^ y < 0, (2.10) 

for some < ;U < 1 y= (least eigenvalue). The latter estimate improves as the initial 

conditions approach the equilibrium point (0, 1, 0); in other words one can pick // closer to 
the eigenvalue 1 — y= by taking e sufficiently small. We will show that these trajectories 
correspond to a singularity of the original metric (j2.ip at a; = 0. 

2.1 Asymptotic analysis at the singularity 

We will be considering solutions of the system (12. Sp . with {W{0),X{0),Y{0)) suffi- 
ciently close to the equilibrium point (0,1,0) and with Y{0),W{0) > 0. (The reflection- 
symmetric trajectories over {Y = 0} and {W = 0} are easily seen to correspond to the 
same metric, while the trajectories with 1^(0) = do not to correspond to Riemannian 
metrics). 

We proceed to derive the asymptotic behavior, as y — t- — oo, of the variables W, X, Y . 
Changing back to x, using (j2.7p . we determine the desired asymptotic behavior of the 
unknown functions in the original system (j2.5p . as a; — )• 0"*". The final estimates will 
confirm that x = is actually a singular point of the metric g. 

Let X{y) = 1 + g{y). Plugging into the equation of W' in (12. Sp we obtain 

W{y) = W{0)exp{y+ [" W{z)g{z){2 + g{z))dz - X C W\z)dz] , 

Jo Jo 

where according to ()2.10p . for y < 0, 



ry I 

/ W{z)g{z){2 + g{z))dz\ < 3e2(2 + V3e)- 
Jo 



e 



2fi 

and 

s/ , 1 _ p3m 



A / W'\z)dz\ < |A| -sVSe^- 
Jo 



e 



Sfj, 



Thus, 



W{y) =Ciey + VF(0)e^ 



exp { / W{z)g{z){2 + g{z))dz - X W^{z)dz} 



exp{-/ W{z)g{z){2 + g{z))dz + X W^{z)dz} 



where Ci = W{0) exp {- J^^W{z)g{z){2 + g{z))dz + X J^^W^{z)dz} > 0. Using dHO]) 
again, we readily estimate the second term as above and obtain 



W{y) = Ciey + 0{e^'^^+^^y) 
Similarly, from the equation of Y' ()2.8p we derive 

Y{y) = Cae^'-^^^ + 0(6^"'^+'"^)^) 



y<o. 



y<o. 



for an appropriate positive {Y{Q) > 0) constant C2. As for X, we already have directly 
from (I21T)]1 the bound 



X = l + g{y) = l + 0{e^y) 



y<o. 



which we can retrieve from the equation of X' by integrating on (—00, y) and using ()2.10p . 
along with the previously derived estimates for VF, Y . Now we wish to derive asymptotics 
for the variables in (12. 5p . According to the transformation (12. 6p . for y < —M (M > 
large), we have 



^ y^^^^^W ^ v/^^(Cie^ + 0(e(2M+%)) ^ y^^^^^^^y ^ ^(e^^^+T^)^) 



y 



C2e^^"^^^ + 0(e^^''+^"^^^) 



ip X 






Co 



V' \/^l^ V^(CieS' + 0(e(2M+i)j/)) ^C7- 



.e-J/ + 0(e('^~i)s') 



1 



V^-i 



uj =nt - — - X^e-y + nO(e'->'-^^y) 



-y + 0(e(/^-%). 



Also, going back to the second equation of (12.51) and dividing both sides by ijP' yields 



± 



{n-l)^ + 



ijj'^ n — 1 ip 



ip^ ip^ ip 



+ -UJ + X 



(n-l)[ 



v^Ci 



+ 0(e(M-i)!/)]^ + 



n — 1 



[i^3p>£le*» + 0(e<*+A''')] = 



+ [7S5r="'+o('"'" 



-i)j/\i r ^^ ^ r-y 






V^ - 1 6-2?/ 



+ 0(e 



(M-2)y^ 



Furthermore, the first equation of ()2.5p gives 



LO = n 



,| - A = -(V^ - 1)^ + nO(e(^-2)j/) + A = -(^^ " 1)^ + ^C 



,(M-2)3/^ 



Having derived asyinptotics, as y — )• — oo, for all the unknown functions appearing in the 
problem, we would like to derive corresponding asymptotics in the independent variable 
X that we started with. For that we recall (|2.7p and normalize so that x — t- 0"^ as y — >• — oo 
to deduce 

x= I Wdy = I Cie^ + 0{e^^^'+^'^y)dy = de^ + ©(e^^^+i^^) (y < 0), 

Hence, it follows 

for y < —M, M > large. Going back to each of the above estimates, we conclude that 
(for X < 1) 

W = x + 0{x^^'-^^), X = l + 0{x''), Y = ^"^^"^ ~ ^^ x^~^ + 0(x^^+^~^), 

a 

^=axv^ +0(x v^ ), -7 = ^-+0(x^-^), (2.11) 

-0 yjn x 



X -(/) n x^ x^ 



1 L 

where a = -^^-^^ — ^-^^ > 0. 

Our analysis verifies in particular that the above trajectories {W,X,Y), emerging 
from the critical point (0, 1,0), correspond to metrics which do not close up smoothly as 
X — )• 0"*". We have proved the following proposition. 

Proposition 2.1 (General existence of singular solitons). For all initial conditions 
(1^(0), X(0), r (0)), y(0) > 0, T^(0) > sufficiently close to the equilibrium point (0, 1, 0), 
the unique solution of I12.8\) exists for all values y G (— oo,0], and corresponds to a Ricci 
soliton metric in the form Ii2.1\) . Along with the associated potential function they solve 
the gradient Ricci soliton equation \2. 2\) [reduced to Ii2.5\) in this context^ and have the 
asymptotic profile 

ip(x) ~ ax^, a > oj{x) ~ as x ^ 0"*". 

X 

These metrics are a priori defined for x € (0,5), for some 5 > small, such that ip^oj 
have a smooth limit, as x ^ 6~ . 

Remark 2.2. In view of the above asymptotic behaviour of ip, as x ^ 0"*", we conclude 
that the gradient Ricci solitons {M'^~^^,g,\7cl), X), n > 1, A G R, that we have constructed 
can be extended to x = as C^ metrics but are singular in C^ norm. 

Remark 2.3. Now we are able to find the limit of the scalar curvature of g, as x ^- 0"*", 
by computing its leading term. Contracting the Ricci tensor of g Ii2.3\) . we obtain 



R{g) = — n— + — 2 (n — 1 — ^-ip — (n — l)^'^ 



1 Tli Tl — 1) X 

2(y/n - 1)^ + -!^ 2-^ - (n - 1)^, as X ^ 0+ (in the sense of (l2lT]) ) 



a^x 
It follows that (for n > 1) 



X n2rj.^ X 



lim R{g) = — oo. 

a:->0+ 
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Remark 2.4. One could also consider the trajectories which emanate from the other 
equilibrium (0,-1,0) of /i2. 8]) (also a source); these in fact correspond to solitons with 
profile 

-il^{x) ~ X ^ u}{x) = , as X ^ 0"*". 

X 

They are in fact defined for all dimensions n + 1 > 2, and in the steady case (A = 0), 
dimension n + 1 = 2, can be explicitly written out as: 

1 2 

i/j{x) = — Lj{x) = , X € (0,+oo). 

X X 

Notice that these metrics are also singular at x = 0, but their evolution under the Ricci 
flow is almost the opposite from the metrics we obtain near the equilibrium at (0,1,0); 
in particular they remain singular for all time. These solitons are however beyond the 
scope of this paper. 

2.2 The half-complete steady solitons 

In the steady case, A = 0, if we consider solutions of (j2.8p with initial point (X(0), y (0)) 
satisfying X^(0) + y^(0) < 1, 1^(0) > and lying close enough to (1,0), the solutions 
of (j2.12p exist for all y € (— oo,+oo). In fact, these trajectories emanating from (1,0) 
translate back to Ricci soliton metrics of the form (12. Ih . which exist (and are smooth) 
for all X G (0, +oo). We will call these Ricci solitons half- complete. This is to distinguish 
them from the general Ricci solitons discussed in the previous subsection, which are only 
known to exist for x G (0, S) for some small 5 > 0. We will focus on the half-complete 
solitons/trajectories in this subsection. 

Fortunately, the known work in the corresponding context of smooth solitons [H 
Chapter 1, §4] can be directly adapted to our setting. In particular, the equation of W 
in ()2.8p becomes redundant and hence (|2.8p reduces to 

(2.12) 




Therefore, in this framework we deal with trajectories in M?. We remark that the 
Bryant soliton corresponds to the unique trajectory emanating from the equilibrium 



point (-^;\/l — :^) and converging (as y — )■ +cxd) to the origin (0,0). On the other 
hand, the source considered in (|2.9p corresponds to the point (1,0). An important tool 
in the analysis of these trajectories is the Lyapunov function 

L = X'^ + Y^, L' = X^{L-l). (2.13) 

Starting (y = 0) with (X(0),y(0)) near (1,0) inside the disc X'^ + Y"^ < 1 with 
y(0) > 0, we easily conclude that the trajectory {X{y),Y{y)) approaches the origin 
(0, 0), as y — )• +00 (at an exponential rate). Whence, it is defined for all y G (— oo, +oo), 
leading to singular gradient steady solitons (Af""''^, y, V(/>), as in Proposition 12. H with 
M""*"^ = (0, +oo) X §"■. One can easily see that the set of all such trajectories fills up the 
domain in the unit disc bounded by the Bryant soliton trajectory (which emanates from 



(-^, -I /l — i)) and the positive X-axis lying inside the unit disc. 
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Remark 2.5. Following JU, Chapter 1, §^y, we derive that the soliton metrics corre- 
sponding to the {X,Y)-orbits above are complete towards x = +00. 

In addition, we get p^, §1.4] that the behaviour oiip{x) and its derivatives, as x — t- +c«, 
is the same with ones corresponding to the Bryant sohton. 

1 1_ _1 ■ _1 _3 ■• _3 

CX2 < ^ < Cx2 cx 2 < -0 < Cx 2 —Cx ^ <ip < —ex 2^ (2.14) 



for some positive constants c, C. The above estimates suffice to derive bounds from 

as X —7- +00, for ah the other variables we will need below. Indeed, we readily estimate 

(for x > M large) 



^ = 0(4,) ^ = 0(4) 1 = 0(1) ij = 0(l 

ip x^ x'^ ip X ip^ x 



^(zi) ^ = 0{-a) V = Oi-) ^2=0{-) -C<c.<-c (2.15) 



Combined with asymptotics ()2.1ip of the trajectories at the singular point x = we have 
the following proposition. 

Proposition 2.2 (Half-complete steady solitons). There exists a 1-parameter family of 
singular steady solitons {M^~^^,g,'V(/)) with M"^^ = (0,+oo) x S" and asymptotic profile 
at X = as in Proposition \2.1\ These solitons correspond to the trajectories (X, Y) of 
\2.1S^) emanating from the equilibrium point (1,0) and lying in the first (open) quadrant 
of the XY -plane inside the unit disk X'^ + Y'^ < 1. Further, the asymptotic behaviour of 
the solitons at x = +00 matches those of the Bryant soliton. 

Remark 2.6. It is worth noting that in dimension five, (i.e., n = 4) the soliton metrics 
and associated diffeomorphisms can in fact be written out explicitly: 

1 6 

il^ix) = a\fx a;(x) = ^, a > 0. (2.16) 

X a^ 

2.3 The evolving soliton metric g{i) and the action of the difFeomor- 
phisms. 



Since the metric g (j2.ip satisfies the gradient Ricci soliton equation ()2.2p , the evolving 
metric g(t) evolves under the Ricci flow 

dtg{t) = -2Ric{g{t)) go = g, 

only via diffeomorphisms. In particular, 

g{t) = {l + 2Xt)-p;{g) te[0,T], 

up to some time T > such that 

e{t) := l + 2Xt >0, 

where 

pt{x,p) = pt{x) pQ = idM 

is the flow generated by the (time dependent) vector fleld 



1 



gradgc, 
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Thus, 

g{t) = e{t) [d{pt{x)f + V2(pi(x))5§"] (2.17) 

We note that since our manifold {M'^'^^,g) is not complete (at x = 0), pt{x) is not 
necessarily defined for all t, but nevertheless it exists locally in time, t G {—£x,£x),x > 0. 
However, as we will see below, by studying the leading behavior of ptix), as x — )■ 0^, for 
x G (0, 5) the flow exists for all future time t E [0, T), for some small T > 0. In particular, 
for the half-complete steady Ricci solitons, which are our main focus here, the flow exists 
for all t > 0. 

Henceforth, we will suppress the sphere coordinates corresponding to different points 
{x,p), {x,q) in M"+^ so that we may consider pt, for each time t, to be a real function 
in X 

pt : (0,(5) -^ (0,+oo), or pt : (0, +oo) -^ (0, +oo) 



Also, abusing notation we identify the time derivative of pt with its coefficient relative 

d_ 

dx 



to the coordinate field ^ component (the others are in fact zero), that is 

dtpt{x) = -77^{^d/dx(t>)pt(x) = -JT^Hptix)) = -^u){pt{x)). (2.18) 



According to our asymptotic analysis (recall ()2.1ip ). 



dtPt{x) = -^^^^ + -^Mptixr-') (2.19) 

e(t) pt[x) e{t) 



for x,t small enough, which also yields 



1 ^/n-l 2 ,/n- 1 

< dtpt[x) < 



e(t) 2pt{x) - ''-'' ' - e{t) pt{x) 
or 

^ , , < 2pt(x)dtpt(x) < ^^ , , — '-. 
e(t) - ^^ ' ^^ ' - e{t) 

Integrating on a small time interval we deduce 

p2(^) + f :/^^dr < pi{x) < pKx) + f ^^^~^^ dT 
^^^ ' Jo 1 + 2Xt - ^t V ; - pov ; j^ ^ ^ 2Ar 

yielding 

/x2 + (V^-l)t, A = 2 /x2+4(V^-l)t 

\ x2 + ^log(l + 2At), A^O -^*^^^-| x2 + 2^1og(l + 2At) ^^'^"^ 

and in particular 

(^-l)t, A = 

/n- 
"2l 



P''tix)>{ '^-1, n^nu, wn (2.21) 
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Remark 2.7. Notice that for all A G M, the lower bounds are strictly positive (when 
n > 1); provided 1 + 2 At > 0. Thus, 

/Ot((0,(5)) C (m(t),+oo) m(t)>0, t>0 

and in particular pt is not surjective. We further note that m[t) is non- decreasing. A 
geometric interpretation of the latter is that the flow pt ''pushes" the domain away from 
the singularity at x = 0. 



We return now to (I2.18P and focus on the singular steady solitons constructed in §2.21 
Besides the estimate (j2.19p which holds close to {x,t) = (0,0), by (|2.15p it follows that 

-C < dtptix) := uj{pt{x)) < -c, (2.22) 

for c, C > and x > M large, t small. Whence, integrating on [0, t] we arrive at 

x-Ct< pt{x) <x-ct, (2.23) 

for an appropriate constant C > 0. 

Further, we can give a complete description of the evolution of the singular Ricci 
solitons (Proposition I2.2p . Indeed, in this case we derive that there is a critical slice 
{^^crit} X S" of the manifold M""*"^ = (0,+oo) x S", which is invariant under the flow 
Pt{)- In particular: 

Uj{x) > 0, X G (0,Xcrit) w(Xcrit) = u{x) < 0, (Xcrit,+Oo) (2.24) 

Thus, for any point x G (0, +00): 

lim pt{x) = XcTit 
and in particular regarding Remark 12. 71 

lim m(t) = Xcrit lim pf ((0, +00)) = [xcrit, +00). 

Since the soliton in this case evolves via the pull-back of the flow /o^ , this implies that 
the limit of the Riemannian manifolds (M^~^^,g{t)), as t — t- +00, is the manifold with 
boundary M^^ = (0, +00) x S" associated with the pseudo-metric Odx^ + ip'^{xcTit)gE,"- 
In order to prove the above picture, it suffices to show that (j2.24p is valid. Then 
the rest follow from a monotonicity argument and the uniqueness of solutions to ODEs. 
From the estimates (I2.1ip . (12.15P we confirm that uj is positive close to x = and is 
negative near -|-oo, hence there exists a point Xcrit where uj{xcTit) = 0. It remains to show 
that this is the only zero of to. We recall at this point a general identity for solutions to 
the gradient Ricci soliton equation (j2.2p (see for instance [6^ Proposition 1.15]). 

Proposition 2.3. Let (M™,(7, V</>) be a gradient Ricci soliton, i.e., a solution of the 
equation 112. S]) . Then the following quantities are constant: 



R + Ag(j) + mX = (tracing) 

i?+|Vg0|2 + 2; 

where R is the scalar curvature of (M™ , g) . 



{ii) i?+|Vg0|2 + 2A0 = Co, 
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Subtracting the two identities of the preceding proposition we obtain 

Ag(/> - |Vg(/>p - 2X(j) + mX = -Co 
Whence, in our context where A = 0, the previous equation amounts to 

4> + t^-<j)^ = -Co, (2.25) 

iP 

Claim: Cq > 0. According to the asymptotics of u; := (j2.1ip . ()2.15p . we easily deduce 
that (j) tends to — c« at both ends of the manifold x — t- 0^, x — )• +oo. This implies that (p 
has a global maximum M, realized at some point x. By (|2.25|) we get Cq > 0. However, 
the constant Cq cannot be zero, otherwise we would have (p = M (by uniqueness of 
ODEs), which of course is not possible. Our claim follows. 

Thus, every critical point of is a strict local maximum. Therefore, (j) can only have 
one critical point, Xcnt = x. ■ 

3 The Stability problem 

Our main goal in this paper is to prove (local in time) well-posedness of the Ricci flow 
for spherically symmetric metrics which are close enough (in certain spaces we construct 
in ^3.2p to the soliton metrics (Propositions 12.11 and 12. 2p we constructed in the previous 
section. We recall below a useful form of the Ricci flow equation for spherically symmetric 
metrics and then proceed to introduce a transformation of our system into new variables 
C,C (13-9I) - These are designed to capture the closeness of the (putative) evolving solution 
under the Ricci flow to the evolution of the background Ricci soliton. The resulting 
system (after a further change of variables, replacing (^ with a new variable r]) involves 
a second order parabolic equation in ^ coupled with a 1st order equation in t], both of 
them having certain singular coefficients. 

The singular coefficients are determined fully by the background evolving soliton 
metric. The precise asymptotics of these coefficients are essential to our further pursuits, 
so we begin by studying those in the following subsection 13. II Next, in §3.2l we set up the 
function spaces in which we will be proving our well-posedness result. In ^4. II we set up 
the iterative procedure by which we will show the existence and uniqueness of solutions 
to the system (j3.12p . which implies the desired existence and uniqueness of solutions to 
the original system (13. 5p and ultimately (13. 2p . 

One final convention: We will be considering perturbations of both the half-complete 
solitons constructed in §2.2l and the general solitons in §2. H and proving the well-posedness 
of the Ricci-flow in both those settingsjj Whenever a distinction needs to be made below, 
the case of perturbations of the half-complete solitons will be denoted by HC, while case 
of the perturbations of the general solitons by G. 

3.1 The Ricci flow for the perturbed metric: A transformed system 

Let ^ be a spherically symmetric metric on M"^^ = (0, S) x S" (-B = 5 > or 

B = +oo), 

a = X^{x)dx'^ + 'ii?{x)gnn, (3.1) 



^Recall that the half-complete solitons are steady, i.e. A = 0. 
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where XjV' ^-re positive smooth functions. Suppose g{t) solves the Ricci flow 

dtm = -2Ric{g{t)) t e [0,T], go = g. (3.2) 

Also, assume that g{t) is of the form 

g{t) = x\x,t)dx^ + iP\x,t)gnn. (3.3) 

We now let s{x,t) be the radial arc-length parameter for the above metric at any given 
time t. In other words, we define 

ds = x{xjt)dx. (3.4) 

Now, expressing ip{-,t) relative to the parameter s (and slightly abusing notation), g{t) 
becomes 

g(t) = ds^ + ijj'^{s,t)gsn. 

As in (j2.3|) . for each t > 0, the Ricci tensor of g{t) is given by 

ti}— 
Ric{g{t)) = -n^ds^ + (n - 1 - V'fe - (n - 1)V'|)5'S"- 

Thus, the Ricci flow ()3.2p . under the above assumptions, reduces to the following coupled 
PDE system: 

f 2xxt = -2{-n^)dsHd„d,) = 2n^x' 
\ 2^Pi,t = -2(n - 1 - i^i^rs - (n - l)V^i) 

or simply 

~ y , ,_^2 x(x,0)=x(2;), ^(x,0)=^(x), tG[0,r]. (3.5) 

V't = V'ss-(n-l)^ 

(Note that the first equation involves the evolution of the radial distance function, while 
the second involves the evolution of the radii of the spheres, at a given radial distance). 

We return now to the metric of the (singular) gradient Ricci soliton (M""*"^, <?, V(/), A) 
(Proposition 12. ip and its corresponding evolving metric (j2.17p 

g(t) = e{t) [d{pt{x)f + V'(Pi(x))<?§n] , 

defined for all t such that e(t) = 1 + 2Xt > 0. Let 

s{x,t) = \/e{t) pt{x), s{x,0) = X ds = \/e(t) dxPt{x)dx, (3.6) 

transforming g{t) into 

g{t) = ds^ + ip'^is,t),gn^ = x^{x,t)dx^ +il?{x,t)gs^, 

where 

X{x,t):=./7{^d,pt{x) i;{x,t):=y^ij{pt{x)). (3.7) 
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Note that tlj{x,0) = ip{x) is the original component of the metric g (j2.ip . Arguing 
similarly to the case of g{t), we find that the evolution of g via 

dtg{t) = -2Ric{g{t)) go = g 

is equivalent to the coupled system 

I ^* ^ "l^^, , , 1-^? X{x, 0) = 1, ^P{x, 0) = ^P{x). (3.8) 

[ Vi =Vss- (n-l)—^ 

We now take a first step towards transforming our system of equations into new 
variables. Let 

C = --l ^ = --1. (3.9) 

X V' 

The above formulas are defined for all x G {0,B),t £ [0, T]. In particular, these variables 
measure (in a refined way) the difference between the unknown functions X) V' &iid the 
background variables X) "0- Note in addition that requiring .^ = at the endpoint x = 
0,t = forces ip to have the same leading order asymptotics at x = as the background 
component tp. 

We next wish to convert (|3.5p into a system of equations for (", ^. We wish to obtain 
an evolution equation in the parameters s and t (s being the arc-length parameter of 
the background evolving Ricci soliton). We are then forced to deal with the discrepancy 
between s,s, which are the arc-length parameters for the evolving metrics g{t),g{t). We 
calculate: 



Ipl 1 X 1 o IMLEU 1 

X XX C + 1 



and hence we write 






Taking time derivatives in (13. 9p and combining (13. 5p . (j3.8p . we derive the following cou- 
pled system for C,C- 



n-, r-, r — n- 



^K + i " '' v(c + i)(c + i) (c + i)(e + i) ^'(c + i) 

Us 

(C + 1)2(^ + 1) 



n 






ij{C + i? {C + i? ' '{C + i?{i + i) i> {C + if {C + if 



It is clear that solving for x, V' iii dSSj) amounts to solving for C,, ^ in the preceding system. 
Now we would like to simplify p.lOp by removing the problematic term nj^r^^rci-Y) from 

17 



the first equation. (Tliis term would not allow us to derive energy estimates of any kind 
for (j3.10p ). In order to do so, we replace the variable ( by 



ry := ^r — —rw- — 1. 



2n 



(e + 1) 



^* = - Mn - 1) ^^ [^^^ - 1] + 2|^^^ii^ 



(3.11) 



The new system of 77, ^ reads 

>s r L 

i-(^ + i)-2 



^Wi+ir!^ e 



^^ 



(e + 1)2"+2 



2n{n - I)- ^"^^ ^' — r/ + 2n(n - 1)^^ ^ 






^■ = '^ + '"-4> L + m + tf-- -«^"]^9^«^'-?TT' <^-^^' 



+ n 



^s 



V'^ '(?? + i)(e + i)^ 



+ 



V' (?? + i)(e + 1)^" (?? + i)(e + 1)^" (r? + i)(c + 1)^"+^ 

1 V's ^s 1 ??s6 



2 V (r/+ 1)2(^ + 1) 



2n-l 



2(?? + i)2(e + i) 



2n' 



r/ 



% 



t=o 



:=eo. 



t=o 



Before attempting to solve the above system, we must first understand its important 
features. It is crucial that we know the exact leading asymptotics of the coefficients in 
(|3.12p . as x,t — ;■ 0+. Recalling that s{x,t) = y^e{t)pt{x) we derive that (for fixed t) 



, , s i l3.7l i i—rT , / , w 13. 6t i—TT , / -5 



i'{s,t) 



1 ^^V^)' 



=(t)V'(- 



U{t) - 



(3.13) 



as long as e(f) = 1 + 2At > 0. Here we recall again that ^(•) is the function in (j2.ip . 
different from tpi-,-) for t > 0. Recall definition (|3.6p and the upper bound (j2.20p of 
p1{x). Going back to (12. lip , for x,t sufficiently small, we deduce 



1 e(t)^ ^ 1 , liL^, lbs 11 . „ 1. 



(3.14) 



^: 



1 1 



+ 0{s' 



/.-2^ 



i^s 



Tp'^ n s"^ Tp ■0^ 



V's IT' 



n 1 



n 



+ 0(s 



M-2N 



Similar estimates hold near (x, t) = (0, 0) for the derivatives of the potential 

uj{s,t) : = 



w(pt(x)) = ^==a;(^==) = ^^ +0(s^ ^) 



dsUi{s,t) 



JM 1 



/6(t) V^W 



(3.15) 



5 . 1 



1 



_ ^ , , (/0t(2;))l = — -w(^^) 

e(t)9,Pi(x)5xLy^ ^^*^ ^^J e{t) ',/7(t)' 



^^ + 0(.^-2). 
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Moreover, as we will see below, we are going to need also estimates for the s-derivatives 
of the coefficients in the first equation of ()3.12p . that is. 



1 , 1 V's en ,e(t)' ^ 1 



6-1- 






U^) =^ - ^ -^ -^4 + 0{s^--) (3.16) 



where b := min{-^,/u}. Further, from (I2.19P and (I3.6p we have 



dts = -^Pt{x) + ^/7{f)dtpt{x) = 4t^ + ^^^^ + 0{s^-^). (3.17) 

It follows that (x, i are sufficiently small) 

2s - s 

and hence 

Vn - 1 < 2s(9iS < 4(Vn - 1), 

Thus, integrating with respect to t we derive 

a;2 + (V^-l)t < s^ <x2 + 4(V^-l)t. (3.18) 

Remark 3.1. According to Remark\2.7\ the lower bound 



inf s^ := b{t) > (^n - l)t (3.19) 

xG(0,5) 



/loMs /or small t > and b{t) is non- decreasing. Hence, we derive 



* 1 /■* 1 / i 

, -dr < -.dT<2. — < +00 (3.20) 

s Jo ^(^-1)t \ y/n-1 



and similarly for the leading term of -r^ 



/■* 1 

/ -^dr < c{t) < +00 (3.21) 

Jo s v^ 



for all n > 1, X G (0,5). However, the time integral on [0,t] of the most singular 
coefficients Ii3.14\ ) of h3.12\) blows up as x — ;■ 0+. Indeed, by (|3.18p it follows that 

/■* 1 1 x2 + 4(Vn-l)t 

lim / ^dT> lim — ^— -log % ^ = +oo. 3.22 

x^o+ Jo s2 x^o+ 4:{^/n - 1) x2 

Now we specialize to the HC case to derive an understanding of the coefficients and 
parameters we are interested in near x = +cxd. First, from (j2.22p . ()2.23p and (j3.6p for 
X > M large and < t < T small we have 

-C<dts<-c c,C>0 (3.23) 
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and 



X 



Ct < s < X — ct, 



(3.24) 



Thus, by (|3.13p . additionally to the above asymptotics (|3.14p . invoking also the estimate 
(j2.15|) we obtain 



1 



v^=^y 7=''V ^=^^^ 



i^i 



1 



^. 



^s 



1 



^ + ("-i)^ = o(^) 



and analogously to ()3.16p 

s4)=o(^) 



<) = o(^) 






^ 

V;2 






Also, corresponding to (j3.15p we have 
-C <uj< -c, 



UJ. 



0(- 



(3.25) 



(3.26) 



(3.27) 



Finally, it will also be useful furtherdown to note a few properties of the less singular 
coefficients in (|3.12p . namely, j^. By (j3.14p and (j3.25p we derive that we can write 



1 

where setting 



A{s,t) 



ds{ 



A{s,t). 






A(t):=P(s,t)||ioo(,), 
see Remark 13. 1[ 



/ A\T)dT = 0{l), 

Jo 



A{s,t) ^^l. 



as t 



(3.28) 



(3.29) 



As explained earlier, the energy spaces we will be dealing with will be defined relative 
to the background arc-length parameter s{x,t) = pt{x). We note some useful formulas 
that will be useful further down in the derivation of our energy estimates. The first issue 
is that due to the coordinate change (j3.6p . the vector fields dg-, dt (the first considered over 
M""*"^, and being the arc- length parameter of the background soliton, while the second 
is defined so that dtx = 0) do not commute. In fact, we find 



d d d 
di'ds ^di 



1 



d 



y^d^Pt{x)dx^ 

A 1 5 1 dtid,pt{x)) d 



+ 



1 



(see again (|3.6 
d d 



e{t)i dxPt{x) dx y/7{^ {dxPt{x)Y dx ' yJ7{i)dxPt{x) dt dx 
A d dxdtptix) d 1 d d 



e{t) ds 
A d 



dxPt{x) ds ^/e(fjdxPt{x)dxdt 



5x[7^w(pt(x))] d Q Q 
e{t) ds dxPt{x) ds ds dt 

x + 6j{pt{x)) d d d 

e{t) ds ds dt 



(using dZlED) 



1 ^ 
e{t) i; 



Pt{x) 



d d d 
ds ds dt 



(from the first equation of (|2.5p ) 
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Thus, using (|2.11|) and (|2.14|) we derive 

[duds] = -^yn- ds 



-2 I r>(^i^-2\ 



[{^ - l)s-2 ^ o{s''-^)]ds, X, t < 1, G 

o\s-'^)ds, {x,t) G (0,+oo) X [0,r], HC 



(3.30) 



Since the Sobolev spaces we wih be considering wih be with respect to the arc-length 
1-forni ds we must also calculate the evolution of this form. 



dtds = dt{\/e{t) dxPt{x)dx) 



A 



W) 



dxPt{x)dx + \J e{t) dtdxPt{x)dx (recall (j3.6 



-— ds + \/e(t) a^ [-^^{Pti^))] dx 



X + uj{pt{x)) 
e(t) 



ds 



1 V 
e(t) ^ 



ds 



Pt{x) 



(using ([2l8]l ) 
(see the first equation of (j2.5p ) 



Employing (j2.1ip and (|2.14p once more we obtain (t > small) 

dfds =^r-: ■ n— ds 

e{t) i, 



/<t) 



[-(V^ - l)s-^ + 0(s'^-^)]a„ X, t « 1, G 



0[s-')d,, 



ix,t) G (0,+oo) X [0,T], HC 



(3.31) 



3.2 The weighted Sobolev spaces and the main result 

As explained the singularities of the coefficients, at (x = 0, i = 0), in the system 
(I3.12p . along with the asymptotic expansions we have derived force us to study well- 
posedness in weighted Sobolev spaces. The weights will be adapted to the singularity at 
X = 0, t = 0. To construct these, we define: 



Definition 3.1. Let 



xo := 



6, G 
1, HC 



(3.32) 



Given a > (to be determined later) and T > small we define 

s^ + at, (x,t) G (0,xo) X [0,T], G 
ip{s,t), (x,t)G [xo,xo + l)x [0,T], HC 
1, (x,t) G [xo + l,+oo) X [0,r], HC 



(3.33) 



where ¥?(-,t) is smooth cut off function interpolating between lim £ (x,i) and 1, for each 

t G [0,T], having bounded derivatives. 

We recall at this point that for the general solitons (G) we consider manifolds in the 
form 



(0,(5) x§", 



(3.34) 
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whereas in the HC case we study the stabihty of the whole half-complete gradient steady 
Ricci solitons presented in ^2.2| see again Propositions 12.11 and I2.2i Further, we shall 
investigate the evolution of the Ricci flow ()3.8p . ()3.5p only for positive time t G [0,r] 
such that 

e(i) := l + 2At > -. (3.35) 

A note is in order here: Since the Sobolev spaces we will be considering are relative to the 
length element ds, it is worth recalling the behaviour of s{x,t) := pt{x) at the endpoints 
X = 0, X = +00 for HC and x = 0, x = (5 for G. We will always consider the values of the 
arc length parameter s at these points to be the endpoints of integration at each time 
t> 0. 

Definition 3.2. Let 

Sminit) ■■= lim s{x,t) so{t) := s{xo,t) G (3.36) 

X-5-0+ 



and 



,{t):= lim s{x,t)^^ +00 HC, (3.37) 



x—^+oo 



for each t > 0. 

Definition 3.3. For any given t G [0,T] we consider the parameter s = pt{x) and let 
H^{s), a > 1Jl\ be the Hilbert space of all functions 

u e H^{s) 

endowed with the norm 

,2 _/■«', , {d^su? 



Mh^^{s) = y ^ + • • • + ;^2^32fc '^S < +00. 

In the special case k = 0, we denote H^{s) by L'^{s). We also let H^q{s) be the closure, 
in H\{s), of all compactly supported smooth functions in (0, +oo) [(0,(5), G] . As usual, 
we define the spaces 

L\Q,T-Hl{s)) L^{0,T;H^^{s)) 

of measurable functions 

u: [0,T]^ H''{s) 

having finite norms 



— / II l|2 J-t II II _ II II 

";H*{s)) •- / IPll^fefs)"^' ll^llL°°{o,T;Hfe(s)) :- esssup||n||^ft(g) 

Jo tPlO.T] 



l7/l|2 



respectively. Often, when the context is clear, we will suppress the parameter s in the 
notation. 



^We are suppressing the parameter t in this notation. 
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0<c<3= o(i)^ <<! ^ ^~' _L ..X „^ (™>1) (3-40) 



We are going to need the following properties of the weight i. By (j3.33p . for fixed 
time t € [0,T], we have 

^ r f, xG (0,xo) 

Qsi = 7l(o,xo) + -^s^illxo.a'o+i) := \ 0(1), X G [xo, xo + 1), HC (3.38) 

[ 0, [xo + l,+oo), HC 

Similarly, for fixed x, employing ()3.17p we derive 

= [^ + j] l{o,xo) + O(l)l[xo,-o+i) (3-39) 

Also, combining (13.190 . (13.240 we obtain the following comparison estimate of s,^ 

"' ' 1 + ^ <[ 1 + 7^' xe(0,xo) 

s2 i ^ ~ I '^' ^^ [a;o,+oo), HC 

for each < t < T. 

Now we can proceed to the well-posedness problem for (j3.12p . First, we introduce the 
function spaces in which we will be constructing the solutions to our problem and then 
state the main theorem. 

Definition 3.4. Assume initially (t = Oj 

Then, two functions rj, ^ defined over M"^^^ x [0, T] are a (strong) solution of /i3.12\) . if 

T] G L^{0, T; Hi) n L\Q, T; H^^^) ^ G L°°(0, T; F^^q) n L\0, T; HI_,^) 

r]teL\0,T;Ll) ^ G i'(0,T; L^^) 

and moreover satisfy the system i3.12\) in the usual sense of test functions. 

We remark here the fact that once we have such a solution to ()3.12p , then we straight- 
forwardly derive that this solution (ry,^) corresponds to a solution of ()3.5p . which in fact 
will be smooth over M"+^ x (0,T], given the parabolicity of the Ricci fiow. 



We define the energy 

Il°°(0,T;H1) "^ II'"IIl2{0,T;H1^i) "^ \\'"\\L°°(0,T;HI) "^ \\''^\\l2{0,T;HI^^) 



£{u,V;T) — \\u\\ i^oc, (nrp rji. + \\u\\r2(nT:H^ ,A + \\'"\\l°°(O.T:HI) + \\'"\\l^(O.T:H^ : ,) (3-41) 



and for brevity let 

^o = l|r/ollii+lieoll^i. (3.42) 

Theorem 3.1. There exist a, a appropriately large andSo,T > sufficiently small, such 
that i3.1^} has a unique solution up to time T 

rieL'^{0,T;Hl)nL\0,T;Hl+,) ^ ^ L'^{0,T;HIo) D L\0,T;HI^,) 

r,t G L2(0, T; L^ ) ^t G ^'(0, T; L^^,) (3.43) 



= m 

t=0 



t=0 
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subject to the Dirichlet boundary condition 







(^and ^ 



for each t G [0,r]. Furthermore the solution satisfies 

£{r,,i-T)<2C£o. 



in case G), 



(3.44) 



(3.45) 



Remark 3.2. In fact it is easy to observe that in the half- complete case the above implies 
that for any T > and a,a > large enough, the initial energy £q can be picked 
sufficiently small to have existence up to time T > 0. On the other hand, the system 
\3.1S^) does not seem to be globally stable (i.e. up to T = ooj in the half- complete case, 
due to the behavior of the coefficients at x = oo. 

4 The Contraction Mapping 

We will prove Theorem 13. II via an iteration scheme, which is essentially a contraction 
mapping argument; see (j4.50p below. We note that throughout the subsequent estimates 
we will use the symbol C to denote a positive constant depending only on n. 



4.1 The iteration scheme and some basic estimates 

We will construct a sequence {??™',C™}„=o s^-tisfying 

im+l j_ ^^"211 (2n\ \f:m\j 



Vt 



"^+1 = 2n(n - i; 



^2 2nr+^+E,t2i7ii^::i_2V'. 



e: 



m+l 



^^ 



Ms,t) ,^ r + 2 






(fm _|_ ly fdm _j_ l\2n+2 






^(■^'^icm ^ +2 ^m+l , V's m+l 



(1 + e 



m\2 



+w 



,1ps 



er' = (— + (n-l 



+ (n - 1) 






-7] 



^+1 _ 2n(r?™ + l)r+^ ET=2{W 



cm\j 



(r/™ + l)(C™ + l) 



271-1 



A{s,t)^^r + 2 



-e 



+ n 



ips 



e: 



m+l 



(^"" + 1 



+ 



i2n-l 



(4.1) 



e: 



m+l 
ss 



^"^ + 1 ijj {ri"^ + 1)(^™ + 1)2" (r/'" + 1)(^"' + 1)2" 



le 



m|2 



iV's 



^, 



,m+l 



1 



'Is Ss 



(t?"* + 1)(^'^ + l)2"+i 2 ?/^ (r/™ + 1)2(^'» + l)2"-i 2 (r/™ + 1)2(^"^ + 1)2" ' 
where we set rj^ = ^^ = and initially 



V 



m+l 



% 



C 



m+l 



i=0 



Co 



m = 0, 1, 



t=o 



Further, ^™+i is required to satisfy the Dirichlet boundary condition 



e 



m+l 



\^min I'^J 



(and C 



m+l 



(so,i) 



0, in case G) 



(4.2) 



(4.3) 



for all ?7i G N, t > 0; recall (j3.36p . (|3.37p . Our main claim from this point onwards will 
be to show that the iterates converge to a solution of (|3.12p . Our method of proof will 
readily imply that any solution of (j3.12p in the sense of Definition 13.41 will be unique. 
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Note that (|4.1|) is linear at each step m + 1, yet the lower-order terms in the RHSs 
associated to the most singular coefficients involve the unknown functions 7^'"+i^^'"+i. 
In view of the asymptotic behavior (j3.14p of these coefficients, as x,t — )• 0"^, and the 
non-integrabihty (in time) of ||-2-||loo(s) (Remark 13. ip . one could not hope to close energy 
estimates for the above system, if one had not set up an iteration of this kind. 

Our irst claim is: 

Proposition 4.1. There exist a, a appropriately large and £o,T > sufficiently small, 
such that 1^. jp has a unique solution at each step m + 1 € {1,2,...}, in the sense of 
Definition \3.4[ In particular: 



V 



GL°^(0,r;/7i)nL2(0,r;/7i+i), r^^ € L'^{0,T;Hlo)nL\0,T;Hl_,^) (4.4) 



Vt 



m+l 



eL\0,T;Ll) 



er^GL^(0,T;L^_i) 



and 



in+l fm+1 



e"'+^GC°°((0,+oo) X [o,r] 



or C°^((0,<5)x[0,r]), for G 



(4.5) 



m G N, r/° = ,^° = 0, subject to ^.2^ , |-^.3| ). Moreover, the sequence {v"^ i ^"^}m=o ^^^ 
uniformly bounded energy 



£{v"',r;T)<2C£o 
for some constant C > 0. 



m = 0,1, . . . 



(4.6) 



We postpone the proof of Proposition 14.11 for the moment to see how we can get 
from this point to the solution of the non-linear system (13.121) . (Actually, we will prove 
Proposition l4.1l in the next subsection by solving more general linear systems of this type). 
It is useful at this point to note the L°° bounds that the functions ry™^™- belonging to 
our weighted Sobolev spaces satisfy. From this point onwards r/"*, ^™', m = 1,2,..., will 
be solutions of (j4.ip (replacing m, + 1 by m) which lie in the energy spaces ()4.4p . 



Lemma 4.1. The iterates rf^,i^,^^ , n G N, obey the following pointwise bounds: 

II '/ Il2 ^ nnn. , i Mr 11? 



^ili-(s,t)<c^c7(fc + i)fo 



Hk 



i.o(,,,)<CC(A; + l)£:o 



(4.7) 



and 



cm 
l^l-|| rcx>/-„^ < C\l C£q\ 






£fc \\L°°(s) — -- V """'uvil^a-i 

for every A: = 0, ... a — 1, a > 1, and a.e. t £ [0,T]. 



tm 
L^+k\m\L^), 



(4.8) 



Proof. For each of r/*", ^"^, ?7i G N, the energy estimate (j4.6p furnishes a point Sv,(t), such 
that 

\^{s.,t)\' + \^—{s.,t)\'<M£o, 

for some M > and a.e. < t <T. Thus, for r/™ and a.e. t we have 



I ' ^ iM2 I '/ /- iM2| 

\-^is,t)\ -\-^{s.,t)\ I 



^^l^l^'^^ 
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thus 



i^ 



{s,t)\'< 



ly-illl n^lll 



^jwif-^irMds 



+1-^(^-^)1' 



(recall (1338]) ) 



< 



C\\'^\\L^{\\'^\\L^+k\\^,\\L^)+2M£, 



Employing ()4.6p once more and similarly for £J^ we obtain (j4.7p . provided /c < q — 1 
We treat ^™ in a different manner for a.e. t G [0,T], tti G N. By assumption ()4.3| 



r{s^in,t) = 0, HC [r (so, t) = 0, G] . (4.9) 

Thus, there exists a point s*(t) G (smin,Smaz] such that S,^{si,,t) = 0. It follows that 

''S cm cm cr, 






s cm 






cm cm cm 

<C|fe||,.(|fe||,. + (fc- 1)11^11,0 






cm, 



< cVC£:o(||^||i2 + (fe - l)|||^||^2), (by (HD) 

for every A; = 0, . . . , a — 1. D 

The iteration (14. ip yields a contraction mapping: 
Definition 4.1. VFe introduce generic notation 

B,D 

to denote rational functions in r/™,^"*, m = 0,1, . . ., satisfying the following conditions: 

• The denomerators of B, D have non-zero constant terms. 

• The constant term in the numerator of B is non-zero, whereas the one in the 
numerator of D vanishes. 

Lemma 4.2. If B, D are functions as above and Sq is sufficiently small, then the follow- 
ing estimates hold: 



0<c< ||fi(s,t)||i^(,) <C 
where k = 0, . . . ,a — 1 and 



.D. 



lL°°(s) 



< CC£q, 



* l|2 I II ^s ii2 



I ^ iM _L II ■^ iM <^ cr'f 



(4.10) 



(4.11) 



for a.e. < t <T and appropriate positive constants c, C; these depend on the coefficients 
of the rational functions B,D. 

We remark that we shall make use of the lower bound of |-B| only to confirm the 
parabolicity condition (see below the second equation of (|4.13p ). 

Proof. The proof is immediate from the pointwise estimates (j4.7p and the energy estimate 

m. n 
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Let 

ar]=r]— r],ak,=c,— c, ?tt, = U, 1,... (4-l^j 

Consider now the two systems (j4.ip corresponding to the steps m + 1 and m. Subtracting 
these two systems, it is straightforward to check that we arrive at a system: 



+ BdCiC + C"^) + IC" Y^e™-B + ^^^^Z)d77'"+^ + 2n(n - l)^dr/™+i 



^^m+i ^ (^ + (^ _ 1)^) [Sdr/^+i + B(ir'+^ + DdC + ^d??"] (4.13) 

+ ^!^Bdr + —BdC+^ + tlB^^idr,"" + dC) + \B\dC^s^^ 

s ip ip 

where we note that each B^D appearing in the preceding system may differ from the 
other. Also, by \B\ in the second equation of ()4.13p above we denote a positive function 
in the B class, having an L°°{s) lower bound c = ^; one can readily check this from (|4.ip . 
()4.6p . provided the initial energy £q is sufficiently small. 

The main goal of this section is to show that the sequence of solutions dr/™"*"^, d^™'"'"^ 
to (|4.13p yield a Cauchy sequence in the energy spaces we have introduced. Specifically: 

Proposition 4.2. There exist appropriately large parameters a, a such that for suffi- 
ciently small £o,T the following estimate holds (for a.e. t G [0,T]j.- 

<|*i?i(r)(||dr?™+i||^, + \\dr^YHl)dr + f B2{T)^\\dr^"^+m,+\\dC"^+m,dT 
+ c(a I A\T)dT + C£o^£{dr]"',d("';T) + CaiC£o)^VT£{dr]'^-\dC^;T), 



where the functions Bi(t),B2{t), t G [0,T] are integrable, and A (t) as defined in /i3.29\) . 

In fact, we will prove part of this estimate in ^4.21 (see (|4.32p ) building up for the 
second half in the end of ^4.31 (|4.44p . We will show in the end of this section that this 
proposition implies our claim, via a standard Gronwall inequality. 

Our goal for this and the next two subsections is to prove Proposition 14.21 Prior to 
embarking on the proof, we outline the main strategy to deal with the singularities of 
the coefficients in the equations (j4.13p . The main issue is that while we start with the 
energies 11 • 11 f/i for df"^, dri"^, we are forced to control their energies in the norms 11 • 11 h-i 
which involve more singular weights; in particular the new weights contain an extra fac- 
tor of -p, for each of the spaces H^ involved. These appear in the LHS of the estimate 
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in Proposition 14.21 All energies || • ||j|^i , || • 11^2 will be called critical energies, and 

the weights j^ or -^ multiplying the functions dS,^ , drj^ , or d^f- , dr]f critical weights^ 
Thus, one theme of our analysis in this section will involve splitting the extra powers of 
1/s among the various terms in our estimates in order to control the resulting critical 
energies as required by Proposition 14.21^1 



We are going to need pointwise estimates for the differences. 
Lemma 4.3. For every m S N and a.e. t G [0,T] the following L°° estimates hold: 



-jtWl^ + (A; + ij II -^^1^x11^2;, 



\L°°{s) 



<Cl 



)> 



(4.14) 



|-^|Iloo(,) < <-UI^^3TllL2+(fc+ij||-^||^2J 



(4.15) 



and 



.drj"^ 2 ^ ^/w^'HT \\2 . n . iMi^^ 



n ^. uUs)<C{\\j^\\h + {k + l)\\^\\h) (4.16) 

+ CT [ \\dC\\l2 dT + CC£oVT [ \\dr~^\\l2 dT, 
Jo °+i Jo "+1 

A; = 0, . . . , a — 1. Further, in the HC case we have the following decay at infinity (for 
a.e. te [0,T]) 



lim dC = 



lim d^^ = 0, 



m = 0,1, . . . 



(4.17) 



Proof. We recall the boundary condition (14. 3p and ()3.38p to obtain 



.dC 



£k 



■is,t)\' 



oJ—:i—\ ds 



^^i^-k%SJ),s 



£k ^fcH 



<C( 



dC' 



\'-f\\h + (k + m^Jh 



) 



for all m G N and a.e. < t < T. We treat d^"^ in a similar manner to ()4.8p . By ()4.3p 
dC^{sminit) = [dC^{sQ,t) = 0, G]. The latter furnishes a point s*(t) e (smm,Smax] 
such that dS,^{Si,,t) = 0, which implies as above 

lf^(M)P<c(ii^lli. + (fc + i)llf^lliO, 

provided fc < a. In order to obtain pointwise estimates for drj'^ we use the equation of 
dryf , analogous to that of d?/™"^ ()4.13p . Setting s = sq := s{xo,t) we have 



_ r^^fo^^ ^ 2n(n - 1)^1 d??'" 
^ So 'V' 



(so,t) 



^i?d^- ^ t±Ddr-' + ^^^^^sdr-'- 



V' 



V^2 



So 



^Bdc + ^sc"^rfC"^ + Bdc~\c~^ + C"^) + ler'T^c""^^ 



(S0,t) 



J"^ = m + l,m,or7n — 1 below. 

The estimate p.40p gives us bounds for -. 
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Note that the coefficients of the preceding equation, at sq := s{xQ,t), are bounded func- 
tions in t. Whence, integrating with respect to time and taking into account the pointwise 
estimates in Lemmas l4.lt l4.2l along with the ones we just proved above we deduce (a > 1, 
a.e. t£ [0,T]) 



.di]"" 



iso,t)\'<CT f\\dC\\l2 dr + CT fwdC-'Wli dr 
Jo "+i Jo °+^ 

-t cm— 1 



+ C\^C£o 



+ CCEqT 



cm— I tm—2 ft 

^s^ II I II '^^■^ II \ ^^ I W^JC 



m-l\\2 
H 



2 dr 



t cm— 2 
II na-1 lli^t 



dt / \\di 



m-l\\2 



1 dr 



< CT \\dCfH2 dT + CCSoVT 

Jo "+^ Jo 



\\dC 



m-l||2 
H 



2 dT. 



(4.18) 



Then integrating in the s-direction yields 



dif' 



is,t)\^< 



f) \^ |2 
Us] 



ds 



+ 



dr]" 



iso,t)\' 



<Ci\\ 



dv. 



^Si\\h + {k+m^\\h)+cTl^\\dn\l2jr 



rT 

+ CC£oVt WdC^Wl^ dT, 
Jo "+^ 



for every m £ N, a.e. t G [0, T] and A; = 0, . . . , a — 1, as required. The decay stated in 
the end of the present lemma follows directly from (j4.4p . D 

Now, let 



^l 



A 



^. 



dF^ := ^DdC + -BdC^ + :!:±BCdr + Bd^iC + C" ) + IC" I^C^ (4.19) 



^ 



i; 



and 



dF^ := {I±t + {n- l)l±){DdC + DdrT) + -BdC + ^B^idrT + d^ 



"^-^ + (n - l)^)(Ddr + Ddii^ + -BdC + ^ 

+ Csidv"" + dc) + Bdcic + er') + ler'i'^c^^"^ + ^r) (4.20) 

+ '^V?B{dr,^ + dC) + BiCdv^ + V?-^dC) + V?-^C-^B{dr]^ + dC) 



In this notation, (I4.13P can be rewritten in the form 



dr\, 



m+l 



€ 

V;2 



V's 



A 



,m+l 



>: 



^^^m+i ^ z±B^^m+L ^ _Ddr^n,+i ^ 2n(n - l)-^(ir?™+^ + dF^ 



A^^m+l 



^C. 



m+l 



,^s 



^ 
^s 



Tps 



^^ + {n- l)^){Bdi]"'+^ + BdC+^) + ^Bdif+^ + \B\di 



^^ 



jrm 



V' 



i; 



V' 



^m+1 

5SS 



(4.21) 



+ ^BdrC+^ + ^^2*" 



We establish some estimates for the functions dF^ , dF^ that we will use in proving 
Proposition 14.21 
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Lemma 4.4. For any function u S L'^{0,T; L\{s)), k = a^a + 1, a>2, a.e. t E [0,r], 
and sufficiently small £o,T the following estimates hold (for k = a + 1 in the first and 
third equation and k = a in the second one): 



u ■ dF^' 
^2° 



■dsdr < Ca 



u 



.dT + C 



+ c( f A\T)dT + C£oVf]\\dC 



Il2 



dr 



m\\2 

L°°{0,r;_ffi)' 



(4.22) 



,^ ^ ^ ^ ^ dsdr < Ca / — ^2^. , ^(iT + C 



U 



jWhdT 



+ -\/C£o 
ei Jo 

ft dm 



m— 1 



fa-1 \\L2 "I" 11 fa-1 11-^^/ 



U 



(ia-1 lli^ 
U 



dT 



pm pm—1 

^ss II I ||?ss II AllJtniii II "■ II j^ 

1 iiij II pa— I wij /ii ^ ii-«^jII/7q,_iii-1j 



+ {CC£o + ei) ||dril|,2^/r + CU / A^(T)dT + C£:o |Mnii-(o,T;Hi 



T 



(4.23) 



and 



f / ^** s ^/' ii^iii^* + f (^^»'^^r ii*'"""?'^./- 



+ jU[ A\T)dT + C£o]£{drj"',dr;T), 



(4.24) 



where C is a constant depending only on n. 



Proof. Recall the estimates on the coefficients (fXTil) . (l3:25]) . (13:28]) . Plugging (i4J9]l in 
the integral below we obtain the following (for a.e. t G [0, T]). 



u ■ dFl^ , 
^—ds 



I2a 



U 



^2 



A 

s 



^ DdC + -BdC + ^BCdC + BdCiC + C~^) + IC^l^dCB 



^ 



(4.25) 
ds 



u ~ d£"^ 

— 1 1 "^11^2 + CCi^oll -^ \\j^2 (using the estimate (j4.10p for the fraction — and B) 



II '^ Il9 ^^ A 



^^^^" l^Wl^ + ll^^lli^ (employ the L°° estimate (jM]), k = Q) 



fia-l 



\ £a+l''L^{Smin,So) 1 1 ^a 1 1 L^ (gg ^ + oo) 



U 

cm cm—l 

?ss II ii,rAmii2 n ?s ||2 

(s)\\ £a-l IIl^ 

2 I r^P^c ^11^'^ l|2 



tm Jtm „, r3 tm cm-1 Jdm 

I n^lnc II ^«« II 11^ l|2 I II ||2 I n^ nc (W ^ss n , nSss II \||"'''S ||2 



I II i|2 I r^. r<c II ''ss II iiJt"i||2 

+ ll^llL2+'^V'-^^oll^^nllL2||de Hi, 



+ C7(^2(^)^^^_j^ 



Sis 



^3S II 
■-Oil ua-l \\L^, 



,2 I /-Y, I r^c m ^ss 

lL2 



+ C\C£ 



^2 (by (lOOll : see * below) 

cm-l ^cm 

_|| I w '^ss II N II "gs ||2 

1 ll-f'^ "1" II ffa-l 11-^^/ II to-1 11^2 



+ C{C£o)-- 



3.. e 



:rllL2i 



.di' 



di 



IL2 



I II ":^s ||2 N 

+ ll7^rTllL2;, 



(utiUzing ([12]) and (gH])) 
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*since i{s,t) = 0(1) when s G (so,+oo); recall Definition 13.11 and the considerations on 
values of s at the endpoints of integration (|3.36p . (j3.37p . Integrating on [0,t], < t <T 
(vT < ^) and utilizing (14.61) once more, we derive (I4.22p . provided CSq < 1. 
We prove now the corresponding estimates for ds{dF]^) 



u-ds{dFY' 

p.a-2 



-ds 



(4.26) 



This time we plug in the RHS of (j4.19p and examine each generated term separately. 
For the first three terms, we recall additionally the estimates on the derivatives of the 
coefficients (I3T6I1 . IJ^TM and (f3:28|) . Then for a.e. t G [0,r] we have 



u 



2a-2 
U 



a. 



%DdC 



,2 



ds 



(^Mh) 



i^l 



€ 



ds{f)Ddr + f,Dsdr + fDde, 



p.a-2 



ds 



— II /iq_i IIl2 + C'Cigoll 2tfa-2 ll-L^ (^y ^^^ pointwise estimate of D (j4.10p . /c = 1) 



's£' 



M 



,de 



L» 



u 



^C" 



I II "• ii2 _j_n\\L2 IK II s ii2 I II "- ii:^ I r'r'p ii_2^_ii-^ 

^^ll^lli2(s_,so) + C||^|li2(,„,+^) (using dSaoD; recah ^M and mM) 
+ CC£oa^\\^\\l2+CC£oa\\-^\\l2 (by dUI]) and the L°° Km ) 



Similarly, utilizing the estimates of B (j4.10p . (|4.1ip we obtain 



u 



a. 



£2a-2 

U 
jp.a-2 

U „2 



A 



-BdC 



ds 



dlSSt) 



ds{-)BdC + -BsdC + -BdCs 
s s s 

dC 



ds (now recall consideration ()3.28p ) 

B.. 



< 11^11^2 + c^^^wii^iii. + ii^iii. + cAHt)\\dn\u,,, . ii^iii. 



« .2 , ^.2..M,rfr 



+ Il77;^lli^ + c^^^(*)ll7^lli^ 



's£' 



< Ca\ 



+ C\\ 



+ CC£oA\t)\\dr\\li +CA^ 
Next, we estimate 



+ <-^ (*)^II^^IIl2 



(by(l33Ql)) 



dC||2 



(invoking the L°° estimate (j4.14p ) 



"" 5 




£2a-2'^ 


d 


l2a-2 



BCdC 



ds 



(SSSb) 



ds{^)Bi^dr + ^BsCdr + jBCsdr + ^^cc 



ds 



s£° 



de 



^ w IM I r<^ r<e w ^^s n n "-"^ \\z 



s^ 



(by (BSD, dnnD) 



M 



+ 1177^1122 + CVC£:, 



'st 






^fQ-l 



IL2 



\\di 



m\\2 



B 



s ||2 



lL°°(s)IUa-l llL2 



n 



n 



I II ll'i I /~'ll^e"i|l'^ II ^ss ||2 I II ^ 11^ I /-'.//^C' II ^ss I 

+ ll^;^^rTllL2 +0||d^ llL-(.)ll-^^rTllL2 + ll^;^^rillL2 + ^.y Gto II ■^^^rTllL2 11-^^^3x11^2 



de: 



s ||2 
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u 



- ^^\\j^\\L-^{Srmn,S0) + C\\j^Zt\\l^(so, + Oo) + ^ V C^OC^Wj^W L^ 






A' 



Il2 



(using (fOnjl ) 



, em 

■^ ■■ Sss II II j(-m||2 



C(C£o)2\\^^\\j^2\\d^"'\\iri (from (HIT]) and the pointwise estimate ^JM) 



cm r^ tm jcm 

I /~(|| ^ss ||2 \\Jcm\\2 I /^. r^c II ^ss n n "''^s ||2 

Employing again the pointwise estimate of i^* ()4.8p . z = m— 1, ?7i, /c = 0, and the properties 
of S ()i30D . (liHD we derive 



u 



jp,a-2 
U 



ds 



m/cm I /-"i— 1"! 



i3der(e+e; 



IC^l'^d^B 



ds 



62Si) 



cmfc'm, I t'Tt— 1^ 



i?.der(c + C"') + BdCsic + c^') + ^^^cic + c. 



^2q-2 

m— 1 



ds 



rz cm tm-L D 

^ n^ np (W ^ss II I ii ^ss II \|| "■ ||2 I ||jcm||2 II -"g ||2 



£1 






iC 



m— 1 



I . II "Sss ||2 I ^ 4 //^c /ll *»«« II I II 'iiss II ^ 



jm— 1 



I II II iiJtTiii l\\ ^ss II I II ^ss II \ I r^. r<c w^^s II II ||2 



iC 



»?i— 1 



m— 1 



I r^w ^iciT^W'i ii^ss iiz I r^. r^c w^^^ w w I|2 i r^. ric 



ICI 



m— 1 



dZ 



s ||2 



to-1 



Il2||-7^3tIIl2 



cm-1 „. rz cm-1 

\^ss II II ||2 I r~<.//^c ll^ss II IIJt»n||2 



I /^^//"T ll^'5'5 II II l|2 I r^. r^c W^ss II ll^tmiU 

V '-'II poL—l 11-^ II ^«— 1 ii-t^ V "II Da—1 11-*^ II ^ ii-t^ 



B 



s ||2 



^(s)ll^„_lllL2 



c 



ei 



iC 



m— 1 



< - VC^^o(||^||.2 + |||^||.2)||— ,||i. + {CC£o + e.)\\dr\\l.^^^ (0 < ei < 1) 



tm—l 



cm cr.o J. „. 

I r^(\\ ^^^ II I ll^s* II All^e"i|i II II 

+ <-ill7^rTllL2 + ||— 3rllL2J||de 11^2^ Jl 7^31 IIl2 



X 



m—1 



iC 



m— 1 



I /^ll'^ss ||2 ||Jtm||2 I r^.lr^c W'^ss II II ^s ||2 , r^( r^c \%\\^ss II 

+ '^l|-^^rrllL2||rfe llj^i +CVC^o||-^^^||L2||;^^rTllL2+C(C^o)2||-^^^||L2||c.^ lU^i, 



(by the L°° estimate (Hl^ ) 

||^^m||2 



applying ()4.14p in the end. (We remark here that the control of the term i?(i^™(^^+^^~^) 
in the above is one of the most delicate that we have to perform here, essentially due 
to the fact that the e nerg ies only involve the first derivative of r/™; this will be used in 
subsection §4.31 below) i^n 

Combining (j4.26b )- ()4.26m . we achieve the following estimate of the term ()4.26p (for 
a.e. t e [o,r]). 



u-dsidFY 



2a~2 



-ds 



(4.27) 



u 



<C(^\\^\\% 

11 «Q. ll_L^ 



+ c\\- 



u 



e: 



c 

m tm—l „. 



— cm cm—1 

<^'?o(||^||l2 + I|-^IIl2] 



-|IL2 



X 



+ c(ii^ii.. + ii^Miidrii^2_^ji^ii.2 + (cc^o+.i)Nriik, 



,2 , /i2/.N_ , ,/^c_^ll ^S II _Ml'^^'' Il2 



+ C{C£,a^ + A\t)a + ^JC£,a\\^U2)\\^\\l2 



1° 



In fact if this term in the equation had been sUghtly more nonhnear, our result would not hold. 
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+ c 



m—1 






e° 



dC 



m 
s ||2 



7° 



-|IL2 



tm m-1 ^ cm c 

I /-i/ll ^ss ||2 I ll^ss ||2 \||Jcm||2 , r<(r'C \% ( W ^^s ii , iiS., 



^.-1...^ ' „^._i„iOiirfniii+c(c£:o)^(ii^ii,. + ii^ii,.)Nni?,i 



rvi— 1 



Thus, integrating on [0,t], t < T < ^, and using the energy estimate (14. 6p . C<?o < 1, we 
obtain the desired bound ()4.23p . 

Lastly, we derive estimates for dF^ (I4.20p in a similar manner to (14.25P . Recall again 
the estimates on the coefficients (|3.14p . (|3.25p . (|3.28p . 



u ■ dF^ 



I2c 






<e\ 



u 



HTIIL2 



' +^C£o 



dC 



-\ — A^ 

e 



dC" 



\g2la-2\\L'^ 

c 



+ 



TIIL2 



(4.^ 



27^32 IIlO (applying (|4.1UP; for f especially) 



st 



•a-l lli^ 



+ 



e 



C7fo" "'" 



?ss II /^ii "'/ ||2 I II "■? ||2 ^ 

(^ r~. cm cm—1 



(by (USD) 



r^ tm (~< r^ tm em- L Jcm 

I {\\^^m\\2 I ||jem||2 \\\ Sss ||2 I ^ . ne (W ^ss II I ||Sss II \||"''^s ||2 

+ j(||dr? ||^«,(,) + \\di ||i^(,))||^^^||i2 + j\J CEo[\\j;^\\l^ + \\j;^\\l^)\\j;;^\\l^ 

m—l 



n rz, cm,—\ 

+ 7VC£(,ll-^^ 



+ ?(ll*''""^ 



L°°(s) 



IL2 



+ 



(lldr? 



m||2 

L°°(s) 



+ ¥i 



m||2 



\ ||?s ||2 

/ ' II Dn-\ lli^ 



s lli°°(-)^ II £«-!■■- e 



C, 

e 



„m-l 
I 'Is ||2 



+ -ms \\l^{s)\\j^^\\l^ + TV '^^ollT^rrllL^ 






<e 



u 



C 



\-^^^UL'^ + jC£oa 



drf' 



+ ^^C£o\ 



em c 

^SS II I II ^, 



a 
m—l 



L2 + 



lis. 



L°°{s)) II ^a_l 

C 

L2 

m—l 



„m /^ r3 em, 

\\\ Is ||2 I "-^ , r^c II Sss ,, 

jll^;;rTllL2 + -V<-^o||^;;zTllL2| 



dr]"" 



s ||2 



^r„2 



„m-l 
I '/s i|2 



Ii2) + j-A^ 



lL2 



L'-(s);il ^Q-1 llL2 



II I II ^s* II \ ^ 

rllL2 + ||— rrllL2j( 



a 



|^^™||2 



Il2 



+ cr 



,(iC™ 



2 , ii^e 
L2 + 



m 
s l|2 



^^^11^2 



+ 



dry! 



^° 



l|2 "l 



+ - 



*-" /II Sss ||2 I (rip 1 



.hXTs-' 



T 






+ (jC£o) 



Wdrrm + Wdrrm 



/o 



m||2 



l^^+ 



dt + cSoVT I \\d^ 

^ Jo 



m-l\\2 
H 



2 dt 



+ ^C£o{ 



dJTs ||2 
/70-1 11^2 



+ 



irferii2 



lL2 



)> 



employing the pointwise estimates (j4.14p . ()4.15p . (|4.16p . Integrating up to time < t < 
T < -^ and using the energy estimate (|4.6p , C£q < 1 , we arrive at (|4.24p . This completes 
the proof the lemma. D 



4.2 The L^ estimates for rf77'"+^ rf^™+^ 

The main energy estimates that will prove Proposition 14.21 are obtained in the present 
subsection and the subsequent one. These are designed to give two Gronwall-type in- 
equalities which together (and combined with Lemma l4.4p imply Proposition 14. 2p . Since 
the reasoning behind the estimates that we are about to derive now and in ^4.31 is quite 
similar, we discuss here the main objective that guides our analysis. As usual, we consider 
the time-derivatives of the L^ norms of drf^'^^ , d,^™'^^ and seek to derive a bound 

^*(N^™"''iii2 + Mr+'iii2) < i?(t)[iidr?™+i|ii2 + iidr+ii2] +F{d7^"\dr,dr,^,...) 

where B{t) is integrable in time. It is precisely at this point that the most singular terma^^l 

^^These are the terms dr;'"+^, d^'"+^ and dr;™+^, d^™+^ multiplied by a coefficient behaving hke ^ and 
^ (respectively) as a;,i — >■ 0^. 



33 



in the system (j4.1|) are very dangerous. A priori, calculating the time derivatives of drj'^'^^ , 
^^m+i -^Q^ijj giyg ijg g^]2 estimatc as above with B{t) behaving like 1/t. It is exactly for this 
reason that we modified the iteration procedure (|4.ip to include the unknown functions 
j^j^m+i^ ^^m+i jj^ specific lower-order terms of the RHSs. As importantly, it is for this 
reason that the weight function £ was introduced. The correct choice of weight function 
is crucial at this step, yielding us an inequality of the form: 

9t(||dr/-+i||i. + WdC'-Xl) < 5(t)[||d77™+i||i. + WdC'-'WhJ +F{drj"^,dr,dvT,.:) 

cm+l\\2 "I 
'S IIl2J, 



+N[\\dv"'+^\\l, +\\dr^^\\l2 +Ner+'|" 



^a + l -^a + 1 



whereas now B{t) is integrable and N is negative, allowing hence the estimate to close. 
This negative term we obtain here will be essential in closing the estimates for Proposi- 
tion [32] in the next subsection. 



We commence with the L^ estimates from dr]"^^^. Taking the time derivative of the 
L^ norm of dr]"^'^^ and using (j3.39p . (j3.3ip we derive (for a.e. <t <T) 



|2 



^ / ^2. ds - acx||-|^||i.(,_,,^) (4.29) 



1^!? ||2 ,rr^\\-l i|2 j-rw-l 



We estimate the last term using the comparison estimate (|3.40p 






since £ = 0(1) in (so,+oo), HC; see (j3.33p and (j3.37p . Using the estimates of the 
coeflttcients H^IM . (132SD, (I3T^ we have 

/ pa ds (plugging in the RHS of d^^^ dOB ) 



£2c 



^5dr+' + ^BdC^^ + -Ddrj"'-^^ + 2n(n - 1)^< 



ds 



n'^??'"+\,9 ^n(^£™+^,9 nf^C^+^,9 C'n(^^"'+^,9 

+ \\^^^-—\\l2 + CC£oA^{t)\\^^^-^\\l2 (by the L°° estimates ofB,D KTOh ) 

^ 7^lllbr lli^(s_,so) + 7ll^^lli2(so,+oo) (employing ([330]), < e < 1) 

^ r^ll d^"'^' l|2 ^ rll ^^"^' l|2 ^ .11 ^^""'' 112 



, ^ii«^ ii2 , / an 0.^1 J 



£«+l lli^(Smm,So) "'" II ^a IIL2(sq^+oo) "■" ^11 OOL 1 1 -L^ 

r+^ II 2 /" d?/"'+^' 



+cc.„^.„ii^iii.+ /^!r::^,. 
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We next consider ^HdC™"*" ||2,2 • Arguing similarly to (I4.29p . we deduce 



1 d 
2dt 



\\d^ 



im+l 



--^■•hil''^ y^'^s-a.""^' "^ 



£a+l ''L^{Sm.in,So) 



(4.30) 



/^ 1 1 ^ 1 1 2 /^~f 1 1 ^ 1 1 2 ^^"f 1 1 '^ 1 1 2 



Again the last term is controlled by (|3.4U|) . so we need only estimate 



d^rn+l^^m+l 



(2a 



ds 



(plugging in the RHS of d^^^ dOTTl ) 



de 



m+l 



(2c 



(% + (n - l)^)(5dr?™+i + 5(ir+^) + —BdC^^ 



+ \B\dC^+^ + 'i^Bdrf^+^ + dFa"^ 



< 



,de 



m+l 



s£" 



+ - 



c ,M 



Il2 

m+l 



+ c 



^"^■'s 
^^m+1 



ds 



(recall (ISia . ([MS])) 



si° 



\h+C\ 



,de 



m+l 



:^m+l 



Il2 



+ 



m+l Jtm+1 



i2 , .,de: 



m+l 



|2 

Il2 



(by (mnil for 5) 



|S|(i^™+MC; 



-ds + 



^ps Bdi'^+^d'qf+^ 



^ »n+l 



(2a 7 V ^^" 



ds + 



di'^+^dK^ 



ds 



+ e 



|2 

Il2 



^2a 

(using dOOjl ) 



fa IIL2(s(,,+oo) 



+ 



-^2^ds 



+ 



I 



m+l Jtm+1 



B\di'^+'di 



ips BdC^+^d'q^+^ 

(2a "■" ' y ^ £2^ 



-(is + 



ds 



We treat the last two terms separately. First, integrating by parts and using the boundary 
condition ()4.3p and ()4.17p we derive 



/ 



\B\di'^+^dif+^ ^ 

(2a ^« 

|^tm+l|2 

l^l '^!. ' ds- 



m+l Jtm+1 



BsdSr+'di 



-ds + 2a 



m+l Jtm+l 



\B\dSJ^+^dC 



(2. 



a+1 



< — W-^ — 



|2 , ,,< 
L2 +'^ 



m+l 



+ e 



der+'„2 , ca^„^^' „2 



'a 
2 Acm+\ 



2 ^c-.de 
^^ + 7 



■m+l 



B. 



( llL°°(s)l 



__f_||2 



jm+1 



+ 






+ 



Ca2^,dC"'+\|2 



e 



IL2 



(see (lim . c= i) 
(recall (l338]) l 



<4ll^ll^^ + ('^ + 7^^° 



+ (-C£:o + 



Ca2 d^ 



m+l 



£Q!+1 lli^{Smin,So) 



+ 



Ii2 (by (lilD) and the L°° estimate (gH])) 



^° 



lL2 



(«>1) 



Similarly, by (j4.3p , ()4.7p , (j4.17p using in addition the estimates of the coefficients (j3.14p , 
mSlM . (I3:25]l . (I3:26]l we obtainEl 



ips BdC'^+'dr]] 



m+lrln^+^ 



v- 



£2" 



-ds 



^^The possibility to control this next term using an integration by parts is essential in order to close 
our estimates for the L^k norms of d^™^^, dr;""^ , without recourse to the higher derivatives. 
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I'4 



us— / — ::^ ds — / — ^t^ ds 



£2q J ^ £2a J ^ pa 



+ 2a I -^ ^ / dsMs (recall ^M. ^M) 



rln"^+^ rlf:^+^ rln^+^ rlf^''^+^ R 

^ W":^ ||2 ■r'll^^^ ||2 ,\\":R ||2 I/^IIZS l|2 II ^s ||2 

- II g£a IIl2 +'-^11 ^^a IIl2 + II g^a Wl^ +'^11 J \\l--{s)\\J^ZT\\l^ 



2 ^ --^ ii "-'/ ||2 ii "-'/ ||2 I 2||^fS ||2 , „2||"'S ||2 

g II g£a "l^ ' II g£a ll-f-^ + "^ II £a+l "L'^{sniin,so) 



+ ^\\ oa IIl2 + tII .oa IIl2 + II „to llL2+a II «c^+i Ili2(, ^ 3 ) + a || .„ ||l2 



< 7^ll^br lli2(s_,so) + 7ll^^lli2(.o,+oo) (employing dllO]), < e < 1) 

+ {Ca + CC^o + «')I|-^IIL(._,.„) (by the L^ estimate dUlD and dMU) 

+ ie + CC£o)\\^\\h + Ca^ll^lli. (a > 1) 

Summing (j4.29p . (|4.30p and taking into account the above estimates we derive that 
for a.e. te [0,T]: 



1 d (^n^+l /7P"^+1 

n\^ m+l||2 I ||Jt'Ti+l||2 \ I fu"''l ||2 , ||"S ||2 \ 

2^ill«^ IIl| + II"'? WlI) + "Cr(|| ^^_^^ \\L^s^i„,so) + II £a+l \\lHs^,„,so)) 



< C7(^ + C7£:oA2(t))(||dr?-+i||2, + H^r+iiO (4.31) 

/^ / *^ \ 1 1 ^ 1 1 2 / 2 \ 1 1 ^ 1 1 2 

+ C(Q + JJII ^^+1 IIl2(s™,„,so) + 7^" + ^Jll £a+l IIl2(s„,„,so) 

+ (4e + -C£:o - ^)\\^\\h + J £2a ds + y ^2. ^^ 

Now we wish to integrate on [0,t], t < T small, use the estimates ()4.22p . (j4.24p in 
Lemma [4.41 and then apply the Gronwall inequality. However, in order to do this we have 
to absorb the terms with critical weights into the LHS. Let e > 0, £^o be sufficiently small 
such that 

lOe + -C£o < T 

e 4 

and a, a be appropriately large satisfying 

1 C \ C 

-a > — -70' > — o- 

4 e 4 e 

Then, integrating on a small time interval [0, t] and invoking (|4.22p . (j4.24p for u equal to 
d7]^~^^ and d^^~^^, e = \ respectively, we deduce 



iWdil'^^'Wh+WdC^^Wl,) (4.32) 



1 
2 



aa f\.drr^, +||^^_||2 Wr + 1 /ll^^^lP dr 

2 / vll £a+l ''L^{smin,so) ~^ '' £a+l ''L'^{smin,so)) ~^ /^ / H £a 11^2 

+ c(a I A\t)dt + C£o]sidr]"\dC;T) + CiC£ofVf£{dr]"'-\dr'^;T) 



< 
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Thus, by Gronwall's inequality and (|3.29|) for t G [0,T], T > sufficiently small, such 
that 

a[ A^{t)dt<C£Q a^T<l (4.33) 

Jo 

we conclude 

esssup(||dr?-+i||i.+||dr+ii2) (4.34) 

ie[o,T] 

< CC£o [£{dC, diT] T) + C£q^/T ■ £{dC^ ,drf~^]T)\ 
and 



aa 







+ \\^^^Tr\\h(sm.r..so))dr (4.35) 



£a + l ''L^(Smin,So)''' £a+l "-^^(«min,So) 



< CCSo [£{dC, drf'-T) + C£Q^/f ■ £{dC^ .dif-^-.T)] . 

4.3 The L2_i estimates for dr/,"^+\ c/^,"^+^ 

Here we seek estimates for the s-derivatives of the terms d^^^^ , drj"^^^ . We remark 
that the estimate (j4.35|) we just derived will be essential in absorbing certain terms that 
appear in this subsection to close the estimates. In particular, it is essential that the 
terms L \\drf^~^'^\\']2 + ||rf'^'"'''^||r2 dt with critical weights appearing in the RHS of 
(j4.4ip have already been controlled in the prior subsection by the energies in the RHS of 
the above estimate (|4.35p . 

Recah (IOqI) . (13:311) and dOO]) to obtain (for a.e. < t < T) 



/■ J„m+1 n J m+l /> I J m+l|2 i i- \jm+l\2 

^ j J2^2 ds -{a- l)a\\ —^ IIl2(,_,,„) (4.36) 



As usual, from (j3.40p 



II "Vs ||2 < r'n-l r ^■^ 11^ I ^11 "Vs ||2 

'' g£a-l \\L'^ — ^^\\ £a \\L^(smin,so) ^ ^ ^ £a-l W L^ (so ,+oo) ^ 

since £ = 0(1), s £ (sq, +oo), HC; see (I3.33p . (I3.37p . In order to estimate the term 

J m+l 

-^dsdvT'-'ds (4.37) 

we plug in the RHS of the equation of drf^^ (j4.2ip and treat each arising term separately. 
For all three of the subsequent terms we use the estimates on the coefficients ()3.14p , ()3.16p , 
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(I22SD, ^M and (I32SD. 



drj. 



m+l 



£2a-2 
P2a-2 



d. 



^. 



2 
s_ 



s R^cm+l 



BdC 



ds 



(B^ZJi) 



—B d£"^~^^ + — 

^2 ^2 ^2 



,dr/: 



m+l 






m+l 



^ II ""'^ ||2 I r*,, ,, 



+ 



dv?^\,2 , ^„^r+\,2 



s£' 



'a-l 11^2 



+ c\V 



ds {B G L~(s) dHUD) 



f_||2 



+ 






di?^' „2 






ndTjg 112 ^^11 '-'''/s ||2 

- ^"^Il £^ \\L'2(s,n^n,so) +^ll £a-l IIl2(so,+oo) 



(employing p.40p ) 






,de 



m+l 



£a+l lli^(Smin,So) 






,der+' 



+ <^ll^^rr IIl2(«o,+oo) (by the L^ estimate gUD and (gUD for S„ C^o < 1) 

Similarly, we obtain 



^2a-2 ^"^ 



-5, 






ds 



g3Zb) 



^^m+l 
^2o-2 



ds{^)BdC^^ + tlBsdC^^ + ^BdCs^^ 






i, 



i, 



ds 



,di 



m+l 



idv. 



m+l 



< 11 ""'^ ||2 I /^i r'^5 ||2 I II "'/^ ||2 I r'll^/J'^+l||2 I 

- II 5£a-l IIl2 +'-11 g^a-1 Wl^ + II g£a-l Wl^ +L'||d^s IIl-^I 



5 



s ||2 
— Ilz,2 



+ e 



,de: 



m+l 



,2 , C'„dr?r+\,2 



-iii^ + 7ii^£^^iii^ 



(recall the L°° estimate (j4.15p l 



<^^||^^s ||2 



'"+^ - C dr]^+^ 



£a lli^(smm,so) e " ^"-1 IIl2(so,+oo) 
I "Ss 

£0 lli^(Smm,S0) """ *^ ll £"-1 Ili^{s0,+Oo) 

Continuing analogously, we derive 



+ Ca\\^\\h,^ . ..^ + C^ll^^^lli2r.„,+oo) + (^ + C7C^o)||^^|li. 



(by dSani), < e < 1) 

dCTs^\ 



drf^ 

pa-2 ^* 
^2a-2 



-Ddr?'^+i + 2n(n - l)%dr/™+i 



^2 

_ 

>2' 



ds 



g37b) 



s s s V 



driT^\,2 , ^^c .2 



+ 2n{n-l)^dri^ 



ds 



„m+l 



s^^ 



dr] 



m+l 



(recall the estimates on D (jiJOj) . (gZU) 
m+l _ _ - „ D, 



s^° 



rlli^ + ll^lli^ + ^^'WMV 



2 , ri /i2/iMi j„m+l||2 



+ 



,t^^! 



m+l 



s£' 



— IIL2 



+ C7^^ 



+ 



,t^^l 



m+l 



dr] 



,m+l 



dflT^ l|2 , II "'/s ||2 ir'll!':!' l|2 I /-ill^'/s ||2 

+ L-ll o«„, 1 llr,2 + O 



s_||2 

L°°(s)ll£a-l IIl2 

^^m+1 



— llL2 -^ II g^a-1 llL2 -^'-'ll^2^Q_illL2 

;^^llL2(s^,„,so)"^'-^ll £a-l IIl2(so,+oo)+C'^ II £a+l "^^(^^^^^^q) (by dMUj)) 

dC+'„2 



j„m+l J„m+1 J„m+1 

^ r<^ll "-^s ||2 I ^11 "^s ||2 _Lr'^2||"^ l|2 



J m+l _ 

C^ll^lli^(.o.+oo)+-C7C7foA^(t)|| ^„ 



2..M|'^^'"+'||2 



i2 + cA 



+ C7C7fo^^(i) 



j„m+l J„m+1 

I "'/ ||2 , II "Vs ||2 

I ^Q. IIl2 "T II £a-l lli^ 



£a-l 11^2 

(utilizing the L°° estimate (|4.16p l 
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+ T [ Wdr^'Wl, dr + CSoVr [ \\drfH2 dT 
Jo °'+^ Jo "+^ 



Combining (14.37b .)-( j4.37b ) with the first equation of (I4.2ip . we obtain the following esti- 
mate for the term ()4.37p . 



^m+l 



^9sdvr'ds 



(4.38) 



c_udv: 



m+l 



< —a 






I (^^|| ^« ||2 I p 2|| "'/ ||2 

' II £« lli^(Smin,So) II £0 + 1 ll-^^(«min,So) 



+ ^-^ll^^lli^(....o) + ^(7 + ^'W)ll¥^lli^ + K^^oA^W + ^)ll^ 



+ c( 



e - ■'■■ t 

La + ll^^^^rrlliO + (^ + '^C'^o;i|-7^3riiL2 



I "s ||2 I II "^s ||2 \ I /, I /^/^c Ml^sss ||2 



^° 



^ ^ drC^^ -WFl) ^^^ 



(2. 



OL-1 



^CCEoA^ify 



T f Wd^^^Wl^ dr + CSoVT [ \\dC\\l2 dT 
Jo "+^ Jo "+^ 



(The term J 



dr,T+'-dsidFl) 



ds will be controlled below by using the estimate (|4.23|) from 



I r2 



Lemma HiU above, putting u = d'rf^'^^). 

We proceed to derive estimates for Hc^C™"*"^ ||^2 • Similarly to (|4.36p . using in addition 

dO]) . (IJTTl) we have 



-jt^ddds + - / '-jtz^dtds 



lil||^A'»+l||2 



/ %T(dsdCr' - n^dC^')ds -(«-!)/ 



< 



1 7tm+l 



£2a-2 



-l^tm+l 



m+l 



L2+CII 



|2 

Il2 



(4.39) 



(a-l)cj||^^^^^^ f.2. +C7(a-1 



dC^'u2 , ^||'^C+\|2 



si° 



riiL2 



There are essentially two terms we must estimate in this case. For both of them, we 
use the estimates on the coefficients (I3.14p . (I3.25p . Recall that by (|3.38p ds£ = 0(1), 
t € [0,T]. Therefore, we start first with the term 



{2a -2) I 



1 r/-m+l 



dCT-^'dCt 



2a-l 



sds (from the RHS of the equation of dC"^^ ([421])) 



m+l 






,t 



>, 



i-^ + in- l)r±)(^Bdr]'^+' + SdC^^) + ^Bdi 



t"^+l^ 



■05 



tm+1 



^ 



+ \B\dCs'^^ + ^Bdrf^+^ + dF^' 



< " „ 

+ C 



2iider+' 



12 ^udv""^\,2 

\2 _^Q^^ I ||2 



7^^llL2 ■^>-^ll^2£a-lllL2 
Jtm+1 ^fm+1 



V' 



ds 



211^"^+' 



V' 



(recall the L°° estimate on B) 



,2 , ^l|rfr+\|2 



a"ll— ^T— 1122 + C||^27^||22 + a 

2 , ,11 ^'^ss ||2 , ^ 2\\dc,s ||2 , 2ii"s, 



2iider+' 



^2£a 



2 
L2 



tm+1 



ea-1 lli^ 



+ e 



— llL2 



H a 



+ C\ 



idv. 



m+l 



+ (2a - 2) / 

c„2l|rfe+^ 



d^^+^dF{^ 



^2a-l 



< —a 



£a ''L^iSmi-n,So) 



sds 

C„2||C"+' 



a "L^ +« II £a llL2 "^'-ll ^^^.-1 Hi^ 

(employ (lOOll . < e < 1; see (l3371l l 



H a 



£a-l IIL2(so,+oo) 
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+ Ca^ 



I £a+l ''L^(smin,so) 



+ CII """'"''■= 



+ Ca\ 






+ Ca' 



,de 



m+l 



£a+l ''L'^{Sm.in,SQ) 



+ c\\ 



d^ 



m+l 



d^'^+l 



^^^^ lli^(Smm,So) ~'~ II £o-l IIL2(so,+Oo) ' ^ £« 1 1 i^ (Smin ,So) 



+ Cct| 



d??: 



,m+l 



J m+l 
+ '-'11 ^o_i IIl2(5q^+oo) 



Jtm+1 /■ 






ods 



We finish the estimates for the -^q_i norm of d^"^^^ by estimating in a similar manner 
the first term in the RHS of ()4.39p . (We remark here that since the term d^^"*"^ must 
necessarily be multiplied by the weight j^^rr in the L^ norms below, the extra powers of s 
arising in the RHS of d^^'^^ K2T\i must all be absorbed in the terms d^"'+\ ^7?""+^ with 
at most one derivative; it is at this point that we make essential use of the estimates on 
the terms d^*""*"^, dri"^~^^ with critical weights that have already been established in §4.2p . 



£2a-2 



-ds 



(plugging in the RHS of the equation of d^f KTHi ) 



£2a-2 



+ 



ip Tp^ Tp 

■Ips 

C..dC^+^ 



+ \B\dCs-^' + ^BdvT^' + dF^ 



ds 



fa-l ''L'2 "T II „2ffa-l ''L^ 



2 , ,|^^s™+'||2 , C.dll 



m+l 



+ £ 



e s^ 
C,dC+\,2 f\B\\di 



Sla^l IIL2 



£2q-2 



'^ ' ds + e -T 



+ — I 
m+l 



,2 , .ndCs^' 



^a~l IlL^^ £ll52£a-lllL2 



+ e 



— IIL2 



+ 



c„dr?; 



m+l 



Sla~l llL2 



^a-1 lli^ 

d^T'dF^ 



fa-2 



ds 



< 4^11 des^+' 112 , c_^2„dr+' 



^^^llL2+7^ 



£a+l ll-f^^(Smm,So) 



+ 



C df 



m+l 



m+l 



^ II £a+l lli^(Smm,S0) 



+ — (J 



+ 



C^drj- 



m+l 



^IIl2(so,+oo) + — '^ 



C,l< 



m+l 



^m+1 



+ 7 II £a-l IIl2(s(,,+oo) + — '^ 



C _udr]s „2 



^a IIl2(so,+oo) 

g_iidei"+\|2 

c„d^f+^,2 



^a ll-L^(Smm,S0) £ " £"-1 ll-f'^{sO, + Oo) 



dC-+ldK 



2 lll^^^lli^ - / fa-2 ds (by dSaOD and (BZUD, c = i, < e < 1) 

Summary: Adding (|4.36p to (j4.39p and combining (j4.38p with the estimates of the 
two terms above we deduce (for a.e. < t < T) 

tm+l||2 



~(ii<i,r'iii;^. + Kr'iiii..) 



+ (a-l)o-( 



I^C+'||2 



+ 



,de+\|2 



£a ''L^{Smin,So) ' H £a ''L'^iSmin,So)J 



(4.40) 



,a 



<C{^ + A\t))i\\dvT+'\\l, +IK 

|rfC+^l2 



+ — (a^ + a) 



+ 



m+l||2 

IIl 

dC+' 



r2 



^a Ili^CSminiSo) II ^i 



his^....))+CC£.A\t)TJ\\dr-''\\l.^Jr 



+ (CCfo + lOe 



iMi^e: 



m+l 



«a-l IIl2 



c 2{ud'n 



m+l 



+ — fJ 



+ 



,de 



m+l 



£a+l lli^(smm,so) II £a+l ll-f'^C 



im )^0j / 



drj-^+'dsidF^) 



ds + (2a - 2) / 

+ C{C£ofA\t)VT [ \\dC\\l2 dr 
Jo "+^ 



m+l J rpm 



d^^+'dFr 



2a-l 



^ds 



02a-2 



-ds 
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We next integrate the above with respect to time to close the estimate, using Lemma [4.41 
and the estimates derived in ^4.2i For this purpose, we need to choose certain parameters 
to be small enough: 

Apply (fi:2i]) with u = -i'^d^"'+^ and let e, i, So be smah such that 



ft I 

CCSo + We + e + CCSqT I A^{T)dT < - 







and a, a appropriately large satisfying 
4(«-l)>- 



4 



1 a2 Q 



in accordance with the considerations after the derivation of (14.3ip . Then, (I4.40p gives 
1 



-Mdv. 



,m+l||2 



I r2 



+ 



(q — 1)(T 



a — l 

,dv. 



(4.41) 



m+l 



+ 



,rfe: 



m+l 



£a ''L^(Smin,So) ' II £a lli^(Smin,S0) 



)dT + - 



1 r,d(z^^ 



4./o " ^°-^ ^^ 



dT 



<C I {a^ + A\t)){\\drj, 





m+l||2 
s r2 



+ \m 



m+l\\2 



)dT 



+ C{T + aC£o / A\T)dT) esssMWdrT^^Wl, + ||(i^-+i||2, ' 

Jo r6[0,T] 



+ Ca' 



tm+l 







+ (2a - 2) 






-)'^^ + i/- 



drj^+^ds{dFY 



£2a~2 



-dsdr 



l2o 



-^isdsdT + C{C£ofVT A\r)dT \\dCfH2 dr 

Jo Jo "+1 



+ c(a f A\T)dT + C8^8{dri'^,dC;T) + C{C£ofVf f 



Ne 



m-l\\2 



ff2 



dr 



We invoke 1) the estimates in the previous subsection (|4.34p . (j4.35p . for T > within 
(I333D, 2) (OHD for u = drif+^ and 3) (H:^ setting u = ^,^(ier+\ e = a io derive for 
a, (T large enough) 



(a — l)cr 






tm+l||2 







+ IMC 



,d< 



'^Ll' 



^(Smin,So) II £a II-L2(s 



.0))^^ + 



t ^cm+l 



,< 



— IIL2 



(4.42) 



dr 



< C 



a^ + A\t) + 



^^o,„ e: 



£1 

tm— 1 



— T r2 + 



,c 



m—1 



0OL-\ 



\l-^> 



(M^: 



m+l||2 



I r2 



MCrii2 )dr 



m+l I 



L? .dT 



ft cm c: 

i n \ i\\ ^s* II I ll^'^s II ^II^C"^ll 

+ '-y^ (ll^^r:TllL2 + ||-^^rrllL2)||de 11^2^^ 
+ {CaCEo^ex)E{drrAC\T)^Ca{CEof^E(dTr-^,dC-^-,T). 
Setting 



ii'i(r) = a2 + ^2(^)^ 



C£:n 



/-m—l 



£1 



tm tr 

X^ZTIIl2-|-|| „^_-^ ||_L2J, 



^2(r) 



tm t' 

'^ss II I II ^ss 

^^3tIIl2 + 



m—l 



(4.43) 



ea-l 



Il2J||^^;— rl|L2, 



7' 
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(lO^ yields 
1 



2^11 Is IIL^_^ 



,._ + mrXij 



(4.44) 



+ 



(q-I)o- 



tm+l 



\ d^s ||2 _|_ II "?s ||2 

I £o lli^(Smin,S0) II ^a lli^(Smin,S0) 



,.o))^^ + 



l ^-^ss ||2 J 



<C/ i7i(r)(||(iC+'f,,2 ,+MC+'f.2 )dr + C/ //2(r)(||ciC+"" 



I r2 






,._^ + ||dC+ii._^)-dr 



/o "-" """ JO 

+ (CaC^o + ei)f (dr?'", dC; T) + Ca{C£ofVT£{dT]"'~\dC^;T). 

It is clear now that we can use some Gronwall type of inequality to get an estimate from 

Lemma 4.5. Xei f : [a,b] ^M be a continuous function which satisfies: 

lf\t) < ^/2 + l^^{T)f{T)dT, t G [a,b], 
where /o G M and ^ nonnegative continuous in [a, b] . Then the estimate 

l\f{t)\<l\M + l^ ^{r)dT, te[a,b] 

holds. 

Applying the preceding lemma to ()4.4ip for 

p2 _ ||j„r?i+l||2 



a.-— I a — 1 



^m+l||2 

'« IIl^ 



^/2 = {CaCSo + ei)£:(<i7?™,dr;T) + Ca{C£ofVT£{dv"'-\dr-';T) 



and 



^(r) = C/7i(r)(||dC+'||i, +||de+i||2, p+CF2(r) 

a—1 a—1 

varying t in [0,T], we derive 



a — 1 



L2 ,. 

a — 1 



(4.45) 



<C H^{r){\\dvT^X^ +^^+'11^ ^dr + C H2{T)dT 
Jo "-1 «-i Jo 

+ {CaCSo + 2ei)^y/ £{dr]"^,dC''^;T) + {CaVf)^C£oy^£{dr]"'-\de''-^\T). 
By ()i:i3]) and (jMD, ([329] we readily obtain 



Hi{T)dT -- 
JO 

H2{T)dT = I 

Jo 



< 



„,2 , a2u\ I V^^° ni ^s^ II I ll^s 

a +A(tjH — {\\j^Zi\\l^ + 



tm— 1 



ea~l 



\Lh 



tm dm— I df™" 

^ss II I II "^sSS II \ II ""^ss II J^ 



£° 



dr < o(l) + 2V2^^^Vt 
(4.46) 



oil ?ss ||2 I oil ?ss 
^IItt — rllr2 "T -^ -rr 



»?i— 1 



-|IL2 



-|IL2 



dr 



c^e: 



ss ||2 



-|IL2 



dr 



< 2V2JC£o ■ ^£{dr^^,di^-T) 
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Thus, from (j4.45p and the standard Gronwall inequahty it fohows that 

esssup{\\dr^T^'\\l, +\\dC^X- f' (4.47) 

Finally, going back to (I4.44p and employing (I4.46p . (I4.47P we derive 

^ + C)exp{o(l) + C-^^^°- 

£1 

+ aVT{C£of£{dr]"'~\dr^^;T) 



(« - 1)^/^ (ll^f^lli^(s_,.o) + \\^^\\hs^..,so))dr (4.48) 

< (o(l) + C^^ + C) exp{o(l) + cVt'^} [{aC£o + ei)£{drr. dST^T) 



and 



/ 



-l^Whdr (4.49) 



< (o(l) + C^^ + C) exp{o(l) + CVf^^} [{aC£o + ei),S(dry'", d^ ; T) 
ei ei 

+ ay/T{C£^f£{d7r-\dC'^;T)\. 

4.4 Proof of Theorem 13.11 

Combining the estimates ([331]), K^ . and (1333) > ([OH]) . ([iliO]) we conclude that 
there exist appropriately small £^^81 such that the relation 



Vf(d77'^+i,d^™+i;r) < K^j£{d'q-^,di-^]T) + K^/£{dr]'^-\d^'^-^;T) (4.50) 

holds, for some < k < i. Consequently, the sequence {??'^,'^'"}^^q is Cauchy with 
respect to the norm ^J£{•, -'^T). It is clear that the latter suffices for convergence of the 
iterates to a solution r/, ^ of (|3.12p . as required in Theorem l3.lt see ([3.4ip and Proposition 
14.11 The energy estimates directly imply the uniqueness of the solution in the function 
spaces under consideration. 

This completes the proof of Theorem 13.11 ■ 

5 The Linear step in the iteration: Proof of Proposition 

In the beginning of ^we took for granted that at each step tti + I, 7tt- = 0,1,..., the 
linear system (|4.ip possessed a (unique) solution with prescribed regularity and energy 
bounds, as formulated in Proposition 14. II We prove these assertions in this section. The 
proof will follow from the study of a general type of such linear systems. The strategy 
for deriving the desired energy estimates in this subsection follows in large part the one 
in ^ Indeed, as we shall see below many derivations are similar and in fact simpler than 
the ones in the previous subsection. Whenever this is the case, we will cite the relevant 
argument in ^ for the sake of brevity. 
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Definition 5.1. We introduce generic notation of functions 

f,geL^{0,T;H\s)) F^ e L\0,T; H\s)), F2 e L\0,T;L\s)) 

satisfying (for a.e. t) 



1 



+ II/IIl-(s) < \\f\\L^{s)+9{s,t) < C 



/. 



s ||2 



for appropriate positive constants c, C , e small, and 
^^ F 



+ 



9s ,,2 



< e <1, 



I 

Jo 



— ||^2dt < +00, i = 1,2 

-ds<Ca\\^ 



(5.1) 



(5.2) 



(5.3) 



u ■ dsFi , ^ ^ II u ii2 



^2i 



:a-2 



£a ''L'^{Smin,So) 



+ G'i(t)||-^|li. + G2(t), 



for the general function u (a.e. < t < T), where Gi{t),G2{t) are positive, t-integrable 
functions. 

Motivated by (14. Ih . we consider the following (more general, as we will see) linear 
system. 



Vt 



V-L^ , V', 



A 



€ 



a + ^/6 - -fr] + 2n(n - l)^r? + F, 



Tp^ Ip 

.Ips 



1P^ 



,€ 



Ips 



Ips 



(t = {^ + {n-l)f,)f-{v + + jfCs+giss + jfVs + F, 



(5.4) 



V 



Vo C 



t=o 



Co, 



e 



t=0 



y^min t^j 



(e 



t > 0. 



(so,*) 



We summarize our goal with the following theorem. 

Tlieorem 5.1. Assume 

Vo ^ H^ ^0 £ H^Q, 

Then, for some small T > there exist a, a sufficiently large such that |5.^[ ) has a unique 
solution up to time T > 0, in the sense of Definition (fg.^D; in particular 



VeL'^iO,T;H'jnL\0,T;Hi^,) ^ e L'^iO,T;Hlo) D L\0,T;HI^,) (5.5) 
Vt€L\0,T;Ll) CtGL\0,T-Ll_,) 



Further, the solution satisfies the energy estimate 



£(ri,i-,T)<C 



io + ET 

Jo 



F, 



« ||2 



T\\i2dt+ / G2{t)dt 




(5.6) 



for some positive constant C = 0(^L Gi(t)df) 



Our first task is to show that the preceding theorem actually implies Proposition l4.1[ 
Throughout the subsequent estimates, we will use the symbol C to denote a (sufficiently 
large) positive constant depending only on n. 
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5.1 Proposition 14. 11 follows from Theorem 15.11 



By the induction hypothesis (see the beginning of the present subsection), we have 
at our disposal the pointwise estimate ()4.7p on the m term of the sequence (f?'",^'"). 
Hence, for sufficiently small initial energy £q, we verify that the system (j4.ip is of the 
form (|5.4p (see Definition 15. ip with 



Fi := 2n(n - 1) 
and 



2 V^2n (2n\ii:m\j-2 
^2 (^m_^i)2n I? I 



^^n^ r + 2 



le: 



m\2 



S (C*" + 1)2 (^™ + l)2"+2 



^2:=-(^ + (n-l)^) 



€.E?"2(?)rp-,,^,2,,„ ..^.„r + 2 



^P 



^2' (-^m_^l)2n- 



— irr + (n-i)-e 



s ^™ + 1 



(5.7) 



(5.8) 



\i 



m\2 



'Is Ss 



{t]"^ + 1)(C"" + l)2"+i 2 (r/™ + 1)2(^™ + 1)2" • 



Indeed, similarly to the derivations of the estimates in Lemma |4.4|, using the bounds on 
the coefficients (l3llD . ^A6\\ . IK25h . ^Mi), dSSS]), the L°° estimates (gZD, (BSD and 
the comparison estimate (I3.40p we can derive the following estimates (for a.e. t G [0,r], 
C£q < 1 and the general function u): 



Fi 



1 l|2 



f)a-l 11^2 
^2 ||2 



to- 



-|lL2 



< c 

< c 



_2^c II? ||2 I „/i2/.M|? ||2 I ^ r^c II Sss II II ?s ||2 

cr C£:nll^r-llr2 +(7A (t)||^::rllr2 + y C£o\\—:^\\l2\\—^\\j^2 



lL2 



lL2 



7° 



^'Cfoll^ll?.+aA2(t)lC"2 






L2 



(5.9) 



+ VG^o| 



c 



to-1 



s II ,'11 Ss ||2 , II Vs ||2 \ 

— IIl2 + ll7^rTllL2J 



lL2l 



and 



u ■ dsFi 



u 



ds< Call — 



+ C\\ 



+ -VG^o||^||l2| 






+ C7(a2Cfo + ^^'(i))fclli2 



(5.10) 



I rii^riC I /l2/'j.\ I . r^c W ^ss n Nn Ss ||2 I ^11 sss ||2 

+ C(aC7fo + vl (i) + V'^^o||-7^nllL2)|l7^rTllL2 +^II7^^IIl2, 



< e < 1. Integrating up to time T and applying (j4.6p . ()3.29p we confirm the first 
equation in (j5.3p : 



T 



F,; 



y" / Ily^lli2rft < C[a2(Cfo)^T + aCe^oiX) + (C7£:o)^yT] , for T sufficiently small. 
JO «^" 

(5.11) 



and 

rT 



I Gi{t)dt <CT + -CSqVt, (5.12) 

Jo £ 

/ G2{t)dt < C[a'^{C8ofT + aCSoo{l) + (C'fo)2\/r] + eC^o (^ > !)• 

JO 
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Thus, the conclusion of Theorem 15.11 yields the desired solution of (|4.ip satisfying the 
assertions in Proposition 14.11 and the energy estimate 



£{r]"'^\C^';T)<C 



<C 



^o + E 



^ Jo 



So + o{l)C£o + eCSo 



< 2C£o, 



(5.13) 



provided e, T are appropriately small. 



5.2 Plan of the proof of Theorem l5.lt A new iteration 

Consider the system 



^* = ^/^ + jf^s + -fv + Mn - l)^V + Fi 



,^. 



>, 



/. , V'. 



Ips 



6 = (^ + (n - 1)^)/ . (r? + + ^^ + —fis + giss + —fris + F, 



ijj ' '^ ' ip'^ ijp ' ip '' ^"^ ' "^^""^ ' ip 



(5.14) 



m 



t=0 



Co, 



t=0 



y^TTiin i'^) 



(e 



(«o,t) 



G) = t > 0, 



where 



ieL'^iO,T;Hlo)riL\0,T-Hl^,] 
is a given function satisfying 



lllf 



L°°(0,T;//l_o) 



+ iieii 



i2(0,T;//2+i) 



<c 



^0 + E 



^\\i,dt + J^ G2{t)dt 



(5.15) 



(5.16) 



for some positive constant C to be determined later. In addition, we assume the estimates 



||^||2 



dt< 



C 



aa^ 



Jo 



F 



I ||2 



-r\\i2dt+ / G2{t)dt 




(5.17) 



and 



e 



s ||2 



^dt < 



C 



^^ + T.[\\-^i\\l-d^ + j[ G2{t)dt 



(5.18) 



r _ 

J^ ll£allL2(._„,,j"''^ (a -1)^7 

We claim that for appropriately large a, o" the preceding system has a unique solution 

r^eL°^{Q,T-Hl)r\L\Q,T-Hl+{) ^ e L^{0,T;HIq) n L^{0,T;HI+^) (5.19) 
r]teL\0,T;Ll) ^t ^ L\0,T;LI_,), 

which satisfies the energy estimates 



£ir],C;T)<C 



Jo 



F, 



^\\i2dt+ / G2{t)dt 





and 



II i ||2 



dt < 



C 



ac^ 



£o + ^£ \\-^\\hdt + £ G2it)dt 



(5.20) 



(5.21) 
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C 

)dt<- -— 

'0'' (a - Ijcr 



^0 + 






£ W^Whdt + £ G2it)dt 



(5.22) 



Observe that if we can prove this, a standard iteration argument (passing to a subse- 
quence, weak limits etc.) yields a solution rj,^ of the original linear problem (j5.4p in the 
same space ()5.19p and satisfying the same estimates as above. This reduces the proof of 
Theorem 15. II to proving our claim above. ■ 



5.3 Apriori estimates for r] 

We proceed to derive estimates for the L^ norm of t]. Utilizing (|3.31|) . (|3.39p we 
derive 



ld_ 
2di 



IvWl^ 



< 



VVt 

£2a 



ds — a 



ds — aa\ 



r?^ 1 

^i^tids + - 

V ||2 



I2c 



■dtds 



+ Ca\ 



V ||2 



(5.23) 



, ri^w ^ ||2 , ^11 V ||2 
First, employing ()3.40p we estimate the last term 



'si- 



V ||2 ^ r<^\\ ^ i|2 I nu ^ ii2 

L2 ^'^<^\\j^^\\ms„- .^>+'- 



o) 



£a IIL2(so,+oo)' 



since i = 0(1), when x £ (xo,+oo); see Definition 13.11 and (|3.36p . ()3.37p . Recall the 
estimates of the coefficients (|3.14p . (j3.25p . (j3.28p . Going back to the first equation of 
(I5.14P and plugging the RHS into the integral below we obtain 



ilVt 



ds 



^. 



V', 



A 



€ 



^fi + ^/6 + -fr] + 2n{n - l)^r/ + F, 



tp' 



V' 



^P^ 



ds 



^ ^11 J^lli^ + II J^lli^ + ll^lli^ (by the pointwise bound on / dS^D) 



+ ^^Wllllli^ + ll7^lli^ 



I V ||2 , (-||ZL||2 



^ n^w ■' I|2 I nw I I|2 I n^ii ^ ii2 

S C-cr|| ,„ , -, \\T2f^ . „.a + (-^llTrllr2,'.„ _L^^ + tvCr 



C 



-4-r\\ ^ i|2 

+ <-^ll£allL2(so,+oo) 



using the comparison estimate (j3.40p once more) 



+ ^ll|lli^+^^Wllllli^ + 



Fi 



Tlli^, 



The last two estimates and ()5.23p yield 

1 '^ II ||2 V ||2 

2dt " "'^ll^a + l lli2(Smin,S0) 

< c(^ + «)ll7^lli2(._,.o) + (^« + ^'W)ll^lli^ + ^^ll7^lli^(.™,„,so) 



Let 



(5.24) 



+ C7| 



6 



^ i|2 -urn'''* l|2 



+ 11^^ IP 
+ IIt^IIlZ' 



-a> C 

4 



-cr > a. 

4 
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By our assumptions on S, (|5.16p . (j5.17p . (|5.18|) . integrating on [0,t], t < T, we deduce 



Mh + 



ao 



^ l|2 



2 / " £a+l "-'^^(■^minjSo) 



dr 



(5.25) 



<l\\m\\li+l {Ca + A\t))Ml,jT + Cal 



I ^ ||2 



dr 



+ C 



'« l|2 



,^1 



£ £ /' 



dr 



provided C is large enough and T small. Thus, by Gronwall's inequality and (j3.29p . 
t G [0, T], we conclude (a large) 



l'?lll,°°(0,r;L2(s)) ^ 



and 



^ l|2 



£a+l Ili^C'^miniSo) 



dr < 



10 



c 



)^^T.[\\^.\\Ur.[a^, 



{T)dT 



aa^ 



Jo 



« ||2 



^||i2(iT+ / G2{T)dT 




(5.26) 



(5.27) 



Now we continue to the L^_i energy estimates of r/^- We are going to need certain 
L°° estimates. Similarly to the derivation of (14. Th . by the energy estimate of £ (I5.16P we 
can obtain the estimates (for a.e. t G [0, T]) 



l|lii«,,)<(k + i)cc 



JO 



^lli2di+ / G2it)dt 

'^ Jo 



and 



l^lli-(.)<^ll^lli2 + ^^c^[^o + E/ llA"'^^' + / ^^^')^' 



(5.28) 



(5.29) 



for every k < a — 1. Moreover, the same argument combined with the above estimate 
dOej) yields 



\j^\\Us)<C\\^\\h + {k + l)CC 



^0+Y. £ \\-^X^d^ + £ G2{t)dt 



(5.30) 



provided k < a — 1. Employing (|3.39p . (|3.30p and ()3.3ip we derive 

1/7 f 71 (j+T) r 71 1 r 71 

■nsdsvt 



-ds + C\\- 



Vs ||2 



(a — l)cr 



'?«l|2 



^2a-2 """ ' "" II g^a-1 11^2 I" ^ H ^a ^^^(s^^^^sq) 
"^s ||2 I /-< II '^^ l|2 I r<|| ''s ||2 



(5.31) 



-I- r'n-H"' 11^ I r'n-ii '^ i|2 I r^ii '« 



s£° 



lL2 



As before, by i^:¥]\i it follows 



I ^s ||2 < /^Q-||!?i||2 I /^ii ^5 ||2 
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In order to estimate the main term we recall additionally the estimates on the derivatives 
of the coefficients ()3.16p . (j3.26p . We will skip the detailed derivations for the following 
estimate, since it is of similar form (only simpler) to the one carefuly claimed for the 
term (|4.37p . In the same spirit we deduce (for a.e. t S [0,T]) 



VsdsVt 

£2, 



a-2 



ds 



(plugging in the RHS of the first equation ()5.14p ) 



f2a-2 



d. 



^/e + ^fis + -h + 2n(n - l)^r/ + Fi 



1 1 ^s ii2 II ^ ii2 ii^ii2 



p-WC ||2 II J'^ ||2 I p| 



s 

Sss II 2 



Js ||2 I II ss II 2 

+ II n„. 1 IIl2 



+ A^ 



si' 



V ||2 I /i2/j.\|| ||2 

tIIl2 +a 1'^JBIIl°°(s)II to3TllL2 



(0 < e < 1) 

Is 11 2 



+ A^ 



II V ii2 1 /• Vs ■ dsFi ^^ 

rlli2 -r II^2^„_iIIl2 + / ^2a-2 "'^ 



r]s ||2 



Now we apply the comparison estimate p.40p . the L°° estimates (|5.28p . (j5.29p . (j5.30p 
and the assumptions (j5.2p . ()5.3p in Definition 15.1! to finally obtain from ()5.3ip 



^ '^ ■■ '12 1 / IN ||%||2 



2 dt ''"^^ "^^-1 ^ ^" " '^^''" F "^'(----o) (5.32) 

< & + c^«)ll^lli^(s_,.o) + [7 + c« + C^^'W + G,{t)] \\Vs\\l.^_^ 



+ C-'\\7^\\Us^,.,so) + [C + CaA\t)] ll.lli. + Ca 



2|| ■? ||2 1 r*!! ^ l|2 



\ £a + l ''L^iSm,in,So) 



-K^n-ll ^'^ l|2 I r'll ^^ ||2 I fi^W ^ss ||2 

+ ^^ll£^llL2(s„,„,.o) + ^ll;^^rTllL2 +'-'^ll7^^rTllL2 
+ G2(t) + [Ca + CA\t)] C So + Yl 



<^ Jo 



Let 



4(a-l)>7 



1 ^ a 

-a > C -. 

4 a- 1 



Then integrating on [0,t], t <T, we derive 

1|| l|2 ,1^ n /■*||^^||2 ^ 

2 a-l Jq e or-l 



(5.33) 



+ C7a' 



^ /^ II ^ lli^(.™„,.o)'^- +{CT + C £ A'{t)dr^ Ml^iO,T;Ll 



(«)) 



? Il2 j_ 1 ^_ / II ?« ||2 



+ '^^' / Il^lli2fs . sn)dr + Ca 



'" 

T 



dT 



+ CTesssnp{\\h\l, + \\^jl,)+Ce [ \\-^X^dT + [ G2{r)dr 

refO.Tl « « Jo ^ Jo 

^0 + E 



+ (CT + C A^{t)dT]C 



^ W^Whdr + r G,{T)dT 
^ JO 
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Recall our assumptions on ^ (|5.16|) . (|5.17|) . (|5.18|) . Utilizing (j3.29p and the derived 
estimates (j5.26p . (j5.27p we conclude that 

1 1 1'^ Tl 

2ll^-lli2_^ + 2(«-l)^y^ \\-f.\\l^{s^.^,s,)dr (5.34) 



< Jl|5xr/o|li2 + f [- + Ca + CA\t) + Gi{t)\hs\\l2 dr 



G2(r)(ir+(o(l) + Ce + -)C 



a' 



^o + E 



T 



Fi 



i ||2 



-\\l2dT+ / G2{T)dT 




Thus, for e, T appropriately small, C and a large, it follows from Gronwall's inequality 
that 



l'?«llL-(0,r;L2_,(s)) 



< 



C_ 

To 






■'^^ l|2 J , 



/ G2(T)d7 

Jo 



and 



T 



£all^'(*™— ^o) - (a_i)(j 



JO 



* ||2 



^||i2dr+ / G2(r)dr 




(5.35) 



. (5.36) 



5.4 Energy estimates for ^: A Galerkin-type argument 

Now that we have solved the first equation of (I5.14p for 77 and obtained the required 
regularity estimates for this solution, we plug it into the second equation of the system 
(|5.14p . We initially seek a weak solution of (|5.14p 

eGL-(0,r;L2(s))nL2(0,r;i?i_,i,o(s)) eiteL\Q,T;H-l,{s)), 

{H'^li being the dual of i:^a+i,o)' satisfying 

^^6,^),,dt =^^ [((^ + (n- 1)||)/ . (, + 0,^),, + (|/6,^^),, 



y 



(5.37) 



^P 



dt 



e 



?o 



e 



t=o 



y^min )!■/ 



(e 



(so,*) 



G)=o tG[o,r], 



for all 

V G L^{0,T;Hl^,{s))nL\0,T;Hl^,is)), 
where by (•, •) r2 we denote the inner product in H^ 

flW2 



(^^1,'V2) 



L2 



I2c 



-ds. 



(5.38) 



(5.39) 



We prove the existence of a (weak) solution to (j5.37p in the energy spaces we claimed. 
(Uniqueness of the strong solution will follow readily from the energy estimates we es- 
tablish along the way). Let {uk{x)}'^^-^^ be an orthonormal basis of L^((0, +00)), which 
is also a basis of Hq ((0, +00)) [(0, 6), G] . Then, for each t G [0, T], 



Wk{s,t) := Tukis- Smin) 



ruk( 



So - Sr 



-)> G 



(5.40) 
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k = 1,2, . . . is an orthonormal basis of L'^{s) [up to uniform rescalling (t fixed), G] and 
a basis of H^ g. We note that by the comparison ()3.40p and (|3.19p 

[ [^dsdt<c[ f\ds<cf -^ds<C I ^dt = o{l) < +00 , (5.41) 

JQ J ^ Jo J S Jq Smin Jo V^ 

as T — )■ O''" . From this we derive that the set of functions 

spa.n{dk{t)wk{s,t) \te[0,T], k = l,2...}, (5.42) 

dk{t) smooth, is dense in L'^{0,T; H^^iq{s)). Recall that by (|3.38p we have dsi = 0(1) 
and hence, similarly to (j5.4ip we verify the asymptotics (as t — t- 0^) 

y^I^^« = 0(^) J ^i^E-d^ = 0{^^) J ^,. ds = 0(-^), (5.43) 

without of course any uniformity in the RHSs with respect to the indices ki,k2 G 
{1,2,...}. Further, by the definition of our basis (|5.4Up we also find 

—^ds = 0{-^^), (5.44) 

since dtl = 0(^"^), as well as dts = 0{s-^ + 1); see (ISlMll . (iSTm . IK23h . 

Now, given a u £ {1,2, . . .}, we construct Galerkin approximations of the solution of 
(j5.37p . which lie in the span of the first v basis elements: 

e := J2 «'^(^)^'= «'^(0) ■■= J ^^^^^0^dx (5.45) 

fc=l 

solving 

(^r,-.)^. = ((^ + (n-i)||)/-(ry+n,-o.^ + (f/er,-.),. 

- {9sC,Wk)^2^ - {gC^dsWk)j^2^ +2a{gCs,Wkj)^2^ (5.46) 

for t G [0, T] and every k = 1, . . . ,1^. 

Proposition 5.1. For each v = 1,2, .. . there exists a unique function ^'^ of the form 
|5.-^5[ j satisfying |5.-^g[ ). 

Proof. Fix a k G {1, . . . ,1^}. We express (j5.46p as a linear ODE system in time and 
estimate its coefficients by plugging into the equation the formula of ^'^ ()5.45p and of the 
basic functions (|5.40p . First, employing (|5.44p we derive 



" 1 

{Ct,Wk)L2 = a'kit) + ^aj{t)OjM{^)- 

7=1 ^ 



Next, utilizing the estimates on the coefficients (j3.14p . (j3.25p . (j5.2p and the above asymp- 
totics (j5.43p we can write 

((^ + (n-l)|i)/-r,-.),.+(|/C,-.),.=i:a,(t)0,,,(-L^^ 
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Recall our assumption Wk E H^^ q and apply (|5.43p once more to estimate 






LI 



5^afc(t)Oj,fc(l) + ^afc(t)Oj,fe(-=), 



where we used the estimates (j5.2p of g and the fact 5^^ = 0(1); see ()3.38p . For the last 
three terms of the RHS of (I5.46|) we set (for a.e. t G [0,T]) 

4(t) := ((^ + (n - 1)^)/ • ^, «^fc)i2 + (^/^^' «^fc)i2 + (i^2, w'OlI 



^ l|2 



^«l|2 



<^ll^lli^+ / 72^^ + ^11^111^+/ TI^^ + CIIt^IIl^ 



i^2 



+ / ^ds 



Observe that d^it) is integrable in time t: Indeed, recall our remark above ()5.4ip and 
our regularity assumption on F^ (j5.3p . The r\ terms are of course t— integrable from the 
estimates we derived in ^5.31 

According to the above, (15.46P reduces to a linear first order ODE system of the form 






afc(i)0i,fc(l)+4(t) 



A; = 1, 



The coefficients of the preceding system are thus singular at t = 0, but they are all 
integrable on [0, T]. This implies local existence and uniqueness of the system and hence 
of Galerkin approximations (|5.45p . (|5.46p of (j5.37p at each step z^s{l,2,...}. D 

We proceed to derive energy estimates for the approximate solutions i^ . 

Proposition 5.2. There exist a, a appropriately large, such that for every v = 1,2, .. . 
the following estimates hold: 



u 



u\\2 



lL-(0,T;L2{s)) 



+ 



,er, 



\hdt 



(5.47) 



C 
< — 

- 10 



^0 + Y. £ \\^\\hdt + £ G2it)dt 



and 



|J1||2 

I £a+l ''L^{Smin,So) 



dt < 



C 



aa^ 



£o + Y. [ \\^i\\h<^^ + [ G2{t)dt 



(5.48) 



for some T > independent v. In addition, 



T \ 2 

{it^v)^2dt 



(5.49) 



C 
< — 
- 10 



^0 + Ef ii#Tiii^^* + f G.(t)<itj^"ii.ii^.^^,, 



^dt, 



for all V = X]fc=i dk{t)wk- 
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Proof. As in (I05D . by ([OT|) and prSUD we have 



I ^ II 2 -L/'^n/ll-^ II 2 I /^^ll ^ 



lirili. < / ^rf«-«<^ll^lli^(s_,so) (5-50) 

I /^ ^ W ^ I|2 I /'^^ll^ ||2 I /'^ii ^ ||2 



Recall that from (j3.40p we estimate the last term as 

||_? ||2 < /^^||_S ||2 _|_r'||S_||2 

using the fact £ = 0(1), x G (0,+cx)); see Definition 13.11 (|3.36|) . Recall the asymptotics 
of the coefficients p.l4p . (j3.25p and the poinwise estimate on f (j5.2p . Since ^ is a linear 
combination (for each fixed t) of the Wfc's, k = 1, . . . , z^, by definition of the inner product 
in L\ and (|5.46p we derive 



< 



r»ii ^ i|2 I r'li ^ i|2 I ^in* i|2 I ^ II ^ ||2 I r, c e^\ 

\ pv fv (J pu 

- 9 II 7^ Ilia + ^11 7^ Ilia + T^^llT^Tilli^ ("^^"§ ^^^ pointwise bounds ^^ oi g) 



_L II ^^ l|2 1/^11^^112 , r>\\ i" ||2 I II -^2 ||2 

< f (^ + «')ll^lli2(._,.o) + f «'llf llia(so,+oo) (utilizing dSaoD, a > 1) 

I ^11 V ||2 -Lr*!!— 1|2 _L r^ii!Zii|2 

I r^ll ^^ ||2 I /'9c- ^^ll'^'^ l|2 I II -^2 ||2 r^ ti' ti'^ 

In order to estimate the remaining term we will need the following L°° bound: As in the 
derivation of (I4.14[) . by the boundary condition we imposed on our space we get 

ii^iii-w<^(ii|ii'+(^+i)ii^in' 

provided < k < a. Now we finish the argument by readily estimating 



|Ss ||2 I 1 ||? ||2 II 9s ||2 

e 



( r, C^ <:'^\ <r c-ll^* l|2 _L 

\9sKs^K )l2^ s ^ii7^iil2 + 7iiyii^°°(s)ii;^^miL2 



<{e + -e 



c.Mie; 



s ||2 



e ■■■«°"^" 



where we made use of the L'^_i estimate on g ()5.2p . Combining the above estimates and 
plugging into (I5.50p we deduce 



2d^ll^1li2+-ll^lli2(._,o) (5-51) 

<-(- + «^)lli^lli2r..,„... + -rili2+(3e + ^e '-'^-'^ 



1^„..„2 , „,,|| e 112 _ 

(J 

£ . '■■J^;^\\L-'(.Sm.n,so) + 7II^''IIl2 -^l-at-r Jt--;|| — 11^2 

I CrrW ^ l|2 -Lrll^'^lP _LnU||2 i rllr) l|2 _l||_^||2 
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Let e, £ be sufficiently small and a, a large satisfying 



C 1 

3e + — e < - 
e 4 



1 c 


1 C 


-a > — 


-cr > —a 


4 £ 


4 e 



Then, integrating on [0, t], t <T, and invoking the estimates on r/ (|5.26p . (|5.27p . (|5.35p . 
()5.36p we conclude 
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(5.52) 
dT 
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1, 



C 



<-,\mk^cjje\\i,dr^c^^ 



^o + T.£\\^\\l^dr + £G2ir)d 
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£o + Y.r\\^\\hdr+ rG2{r)dr 



T 



F2 



+ 1^ Wj^rWhdr, 



'0 '^ JO 

where ^g •= ^'^{x,0)- Since {wk}kez is an orthonormal basis of L^, from (|5.45p it follows 

ii^oiii^ < ii^oiii.. 

Thus, by Gronwall's inequality, t G [0,T], we have (T > small, C large) 

rT cu 



lie 



uii2 



lL-(0,T;L2(^)) 



+ 










lL2' 



(5.53) 



<C7 



C 



ao" 



-^o + E/ W^Wh^^ + l G2{t)dt 



For a sufficiently large, together with (|5.52p . we achieve the required estimates (j5.47p . 
(j5.48p . Finally, setting v = Yl'k=i'^k{t)wk and going back to (j5.46p we apply the same 
arguments as above to obtain the estimate (for a.e. t E [0,T]) 



(er,-).,<c(ii-^ib + iiSM 



+ « lll^lli^ + 



I ' II ill ^ II I II ^s II (K t:.a\ 

l^2^^3TllL2 + ll^2^^3TllL2 + ll77^nllL2 (5-54) 



's^° 



s(.' 



ris F2 

rllL2 + ||— 3tIIl2 



(5.55) 



Integrate now on [0, T], apply C-S and take squares both sides. Employing first the 
comparison estimate (I3.40p . then the estimates on rj (15.26p . (15.27|) . (I5.35p . (I5.36P and the 
just derived estimates of ^'^ (|5.53p . ()5.48p we arrive at ()5.49p . for a, a large enough. D 

The estimates in Proposition 15. 21 suffice to pass to a subsequence (applying a diagonal 
argument due to (|5.49p ) and take weak limits in (j5.46p . after integrating on [0,T], to 
obtain a weak solution ^ of (|5.14|) , in the sense of (I5.37|) , and in the desired spaces ()5.56p 

^eL^{0,T;Ll{s))nL^{0,T;Hl^^^o{s)) f^teL\0,T;H-l^{s)), (5.56) 
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satisfying the energy estimates 

cT 



and 



m 



L°°(0,T;L2(s)) 



+ 
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< — 
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dt 



^o+E 



W^Whdt + l^ G2{t)dt 
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l^a+1 lli^(Smin,S0) 



dt< 



C 



aa^ 



■^o + E/ \\^\\hdt + l^ G2{t)dt 



(5.57) 



(5.58) 



5.5 Improved regularity and energy estimates for ^ 



We now show that ^ is a strong solution of (|5.14p . Consider a positive time Q <tQ <T. 
Looking at the second equation of ()5.14p for t G [t^^T], we notice that those coefficients 
which involve V' and its derivatives are smooth and bounded, while f,g€ L°° (0, T; H^ (s) ) 
by consideration ()5.ip . Moreover, from §5.3l we have t] S L°°(0, T; H^) and by assumption 
Fi S L^(0, T; L^(s)), i = 1,2. Hence, by standard theory of parabolic equations [131 §7.1 
Theorem 5], the solution ^ of (I5.14p that we established above (I5.56p . having "initial 
data" C('5;io) £ H^{s) (for a.e. < to < T), attains interior regularity 



CeL'^{to,T;H^{s))nL^{to,T;H^{s)) 



it^L\t^,T-L\s)) 



Since to £ (0, T) is arbitrary, the following derivations make sense. Recall that for fixed 
t > 0, the weight l"^ is bounded below; see Definition 13.11 
We proceed to derive the regularity 



UL°°{^.T-Hl^)rM\Q,T-Hl^{) 



^teL\0,T;Ll_,] 



and the remaining energy estimates we claimed in §5.21 for our solution ^ (I5.56P of (I5.14p . 
Similarly to ([09]l . utilizing the boundary condition in (f5l^ and ([OO]) . (fOTT) . (l330]l . 
^:M^ we derive (for a.e. t G [0,r]) 
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(5.59) 



(a — l)o"| 



£a ''L^{Smin,So) 



+ C(a-l)|||||i.+C||-^ 



l|2 



As usual, we replace the weight s with £ employing the comparison estimate ()3.40p 

II ^s ||2 ^r'^ll^*l|2 I r<\\ ss ||2 



There are two terms we need estimate. For what follows we recall the estimates on the 
coefficients (|3.14p . (j3.25p . Also, note that by (j3.38p we get the crude (but sufficient) 
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estimate ds^- = 0(1). 
• (2a — 2) / p2a-i ^^^^ (plugging in the RHS of the second equation of (j5.14p l 
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(using the -L°° estimate (15.2 
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(0 < e < 1) 



Next, we estimate the term 
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(substituting from (j5.14p l 
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(employing (15. 2p again for /, g) 
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(applying ([330])) 
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Plugging the above estimates into (|5.59p we deduce 
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Let e be small and a, a large such that 
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3e < - 
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4 e a- 1 



Then, integrating on [0, t], t < T, and invoking the derived estimates in the previous 
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subsections for 77 (1S2SD, (IS12ZD, dSlMD, dOHjl and ^ (I071) . (I05|) we obtain 
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Thus, Gronwall's inequality, t G [0,T], we finally conclude (for T > small, a large) 
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This concludes our proof of the linear step. 



(5.62) 
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